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ABSTRACT 

Radar  target  identification  using  electromagnetic  scat- 
tering in  the  Mie  or  resonant  region  is  discussed. 
Frequencies  being  used  are  those  corresponding  to  wavel- 
engths from  half  the  target  size  to  one-tenth  the  target 
dimensions.  The  frequency  domain  calibration  of  the  free 
field  radar  range  at  the  Naval  Postgraduate  School  is 
performed  and  validated  on  several  perfectly  conducting 
metallic  spheres.  Measurements  of  the  radar  cross  section 
and  the  phase  of  the  return  on  several  finite  length,  thin 
walled,  right  circular  cylinders  were  done. 
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I.  INTRODUCTION 

A.   TARGET  IDENTIFICATION 

In  most  radar  applications  so  far,  the  only  properties 
of  the  target  that  are  measured  are  its  location  in  range 
and  angle.  Such  radars  are  sometimes  called  blob  detectors 
since  they  recognize  targets  only  as  "blobs"  located  some- 
where in  the  space.  In  the  era  of  modern  electronic  warfare 
it  is  desirable,  however,  to  extract  more  information  about 
the  target.  Radar  should  be  able  to  recognize  one  type  of 
target  from  another;  that  is, to  determine  that  the  target  is 
a  Boeing  747  aircraft  and  not  a  DC-10,  or  that  the  partic- 
ular ship  is  a  warship  and  not  a  freighter.  This  capability 
is  called  target  identification  or  classification. 

Generally,  radar  target  identification  requires  a  larger 
signal  to  noise  ratio  then  normally  needed  for  detection. 
Hence,  it  is  to  be  expected  that  the  range  at  which  target 
identification  can  be  made  would  be  less  than  the  range  at 
which  the  target  can  be  first  detected. 

One  of  the  approaches  that  might  be  used  for  target 
identification  is  that  of  electromagnetic  back  scattering. 

The  frequencies  that  are  suggested  as  being  important 
are  those  corresponding  to  wavelengths  from  half  the  target 
size  to  wavelengths  about  one-tenth  the  target  dimensions 
The  use  of  such  frequencies  is  said  to  provide  overall 
dimensions,  approximate  shape,  and  material  composition  of 
the  target  [Ref .  1] . 

It  is  particularly  useful  to  investigate  the  resonances 
excited  by  electromagnetic  fields  of  different  frequencies 
incident  on  the  target.  At  resonant  frequencies  the  scat- 
tered   fields   are    stronger    compared   to    non-resonant 
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situation.  The  back  scattering  cross  sections  are  larger 
and  provide  easier  target  detection  and  identification. 

Associated  with  resonant  frequencies  are  the  natural 
modes.  A  natural  mode  is  defined  as  a  possible  "source- 
free"  distribution  of  surface  current  on  the  scatterer 
according  to  Stratton  [Ref.  6].  When  the  target  is  irradi- 
ated with  an  incident  electromagnetic  wave,  surface  current 
is  excited  on  it.  This  surface  current  radiates  and  gener- 
ates the  scattered  field.  A  resonance  occurs  when  one  or 
certain  combination  of  the  natural  modes  responds  very 
strongly  to  a  particular  type  of  excitation.  The  natural 
modes  are  defined  only  by  the  geometrical  structure  and  the 
composition  of  the  target,  hence  the  resonant  frequencies 
will  be  only  functions  of  the  target's  physical  configura- 
tion. Note  that  the  resonant  frequencies  refered  to  in  this 
thesis  are  real  constants,  in  contrast  to  the  "natural  reso- 
nance" frequencies  of  the  Singularity  Expansion  Method 
[Ref.  2]. 

Investigation  of  target  back  scattering  over  a  broad 
frequency  band  appears  to  be  a  promising  method  to  achieve 
target  identification.  If  the  phase  shift  as  well  as  the 
amplitude  of  the  scattered  field  are  measured,  fewer 
frequencies  might  be  utilized  than  when  amplitude  alone  is 
obtained  [Ref.  3] . 

The  polarization  of  the  radar  back  scattered  energy 
depends  on  the  target  properties  and  differs,  in  general, 
with  the  polarization  of  the  energy  incident  on  the  target. 
This  property  can  be  used  as  a  possible  basis  for  discrimi- 
nating one  target  from  another  [Ref.  4].  For  example,  a 
thin  straight  wire  can  be  readily  distinguished  from  a  homo- 
geneous sphere  by  observing  the  variation  of  the  back  scat- 
tered field  amplitude  as  the  polarization  is  rotated.  The 
field  amplitude  from  the  sphere  will  be  unmodulated,  and  the 
back   scattered  field   amplitude   from   the  wire   will   vary 
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between  a  maximum  and  minimum  at  twice  the  rate  at  which  the 
polarization  is  rotated.  By  transmitting  two  orthogonal 
polarizations  and  measuring  the  amplitude  and  phase  of  the 
received  back  scattered  field  on  each  polarization,  as  well 
as  the  cross-polarization  component,  additional  means  of 
target  identification  can  be  provided.  It  is  worthwhile  to 
indicate  that  non-cooperative  target  recognition  through 
various  methods  has  been  an  ongoing  topic  of  intensive 
research  for  at  least  two  decades . 

B.   THESIS  OBJECTIVES 

This  work  should  serve  as  the  basis  of  what  is  a  long 
term  ongoing  project  at  the  Naval  Postgraduate  School.  It  is 
the  theoretical  and  experimental  studies  of  the  scattering 
characteristics  of  targets  and  the  development  of  radar 
target  identification  schemes. 

The  main  goal  of  this  work  is  to  calibrate  the  contin- 
uous wave  (CW)  scattering  range.  This  is  the  first  step 
toward  studing  target  scattering  characteristics  for  target 
identification.  Calibration  has  been  performed  by  using  an 
assortment  of  perfectly  conducting  metallic  spheres  of  known 
cross  sections.  The  theoretical  solution  has  been  tested 
against  experimental  results  for  the  spheres  of  different 
diameters  and  two  orthogonal  polarizations,  horizontal  and 
vertical.  After  calibration,  the  next  step  was  a  series  of 
measurements  on  the  right  circular  cylinders  having  very 
thin  walls  and  finite  lengths  made  of  brass  and  stainless 
steel.  They  come  in  three  diameters  and  several  different 
length  to  diameter  ratios.  Measurements  of  magnitudes  and 
phases  of  back  scattering  fields  from  the  tubular  cylinders 
were  carried  out  and  corresponding  back  scattering  cross 
sections  are  obtained.  Those  results  should  serve  as  the 
basis   for  future   works  in   the  ongoing   project.   A   right 
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circular  cylinder  is  chosen  because  of  its  resemblance  to  a 
missile  body.  Progressively  more  complicated  models  evolving 
from  a  tubular  cylinder  to  a  missile  will  be  constructed  and 
studied.  Each  model  which  represents  a  step  forward  in 
increased  complexity  from  a  simpler  one  will  be  considered 
as  a  perturbation  to  the  physical  structure  of  the  simpler 
model.  By  comparing  the  scattering  data  of  the  models  to 
those  of  tubular  cylinders,  the  effects  of  the  successive 
perturbations  to  the  physical  structure  on  the  back  scat- 
tering cross  section  and  phase  shift  can  be  observed. 
Explanations  of  these  effects  is  currently  being  attempted. 
Application  of  this  technique  toward  target  identification 
and  radar  cross  section  reduction  will  be  explored. 

In  chapter  II  the  general  approach  to  the  solution  of 
the  problem  of  electromagnetic  scattering  is  given.  In 
particular,  the  exact  solution  for  the  scattering  by  a 
perfectly  conducting  sphere  is  presented  along  with  the 
considerations  on  the  back  scattering  cross  section. 

Chapter  III  describes  experimental  setup  in  the  scat- 
tering laboratory  of  the  Naval  Postgraduate  School  and  the 
measurement  procedure,  as  well  as  the  validation  of  measure- 
ment data. 

Computer  programs  for  the  evaluation  of  the  data  for  a 
sphere,  system  calibration  and  target  data  acquisition  are 
described  in  the  chapter  IV. 

Chapter  V  contains  a  summary  and  conclusions.  It  also 
discusses  results  and  problems  encountered  during  this  work. 
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II.  ELECTROMAGNETIC  SCATTERING  PROBLEMS  AND  SOLUTIONS 

A.   GENERAL  APPROACH  TO  THE  PROBLEM 

A  typical  problem  in  electromagnetic  scattering  involves 
three  principal  elements;  these  are: 

1.  A  transmitting  system  consisting  of  an  antenna  or 
radiator  that  is  excited  by  a  generator  of  an  electromag- 
netic disturbance  of  known  behavior  in  time; 

2.  An  obstacle  of  arbitrary  shape,  size,  orientation, 
and  material  that  is  fixed  or  is  moving  in  any  direction  and 
manner  at  an  unrestricted  distance  from  the  source; 

3.  Receiving  antenna  with  associated  apparatus  for 
determining  amplitude,  phase,  and  the  polarization  of  the 
electromagnetic  field  at  any  point  in  space,  including  the 
surface  (or  surfaces)  of  the  obstacle. 

The  oscillating  charges  and  currents  in  the  driving  or 
primary  antenna  constitute  the  source.  The  electromagnetic 
field  maintainedby  them  is  the  primary  or  incident  field.  It 
is  denoted  by  a  superscript  "i":  thus,  E  ,  H  .  The  parasitic 
element  under  study  is  known  as  obstacle,  radiator,  scat- 
terer  or  target,  depending  partially  on  its  shape  and  func- 
tion. The  current  in  the  source  induces  time-varying 
distributions  of  oscillating  charges  and  currents  in  the 
scatterer.  These,  in  turn,  maintain  an  electromagnetic  field 
that  is  known   as  a  scattered,   or  reradiated   field.   It  is 

"^S    ~*s 

denoted  by  a  superscript  s  :  thus,  E  ,  H  •  The  total  field 
E  ,  H  ,  is  the  vector  sum  of  the  incident  and  the  scattered 
fields.  It  is  a  resultant  field  at  any  point  when  the 
obstacle   is   present   and  the   source   is   active.     Thus, 

Ex=  e  +  E,  H  =  H  +  Hs  ,  Figure  2.1. 
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In  general,  the  incident  field  is  not  independent  of  the 
scattered  field  in  the  sense  that  the  currents  and  the 
charges  in  the  source  are  not  the  same  when  the  obstacle  is 
present  as  when  it  is  removed  to  infinity.  The  incident 
field   does   depend  upon   the   presence   and  nature   of   the 
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Figure  2.1    Configuration  of  a  Scattering  Problem 

obstacle,  at  an  arbitrary  point  in  space.  In  order  to 
simplify  the  theoretical  problem,  it  is  usually  assumed  that 
the  source  is  so  loosely  coupled  to  the  obstacle,  usually  by 
a  sufficiently  large  separation,  that  the  currents  and 
charges  in  the  source  may  be  regarded  as  independent  of 
those  in  the  obstacle.  In  this  case,  the  incident  field  at 
any  point  may  be  treated  as  a  quantity  that  is  independent 
of  the  scattered  field;  it  may,  therefore,  be  determined  in 
the  absence  of  the  obstacle.  Alternatively,  the  scattered 
field  may  be  obtained  by  subtracting  the  incident  field  in 
the  absence  of  the  obstacle  from  the  total  field  at  the  same 
point  when  the  obstacle  is  present . 
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The  fundamental  problem  in  scattering  is  the  determina- 
tion of  the  scattered  field  amplitude,  phase  and  polariza- 
tion. The  scattering  of  EM  wave  from  an  obstacle  is  a 
boundary  value  problem.  It  is  often  hard,  if  not  impos- 
sible, to  obtain  the  exact  solution,  unless  the  shape  of  the 
obstacle  coincides  with  one,  of  the  limited,  simple  coordi- 
nate systems . 

To  solve  the  problem,  it  is  necessary,  as  usual,  to 
start  from  Maxwell  equations.  The  first  step  is  to  determine 
the  scattered  EM  field  in  the  surrounding  space  due  to  the 
currents  on  the  surface  of  the  scatterer.  The  boundary 
conditions  on  the  surface  are  then  applied  to  obtain  a  set 
of  coupled  equations,  containing  all  components  of  the 
surface  current.  Only  recently  has  become  practical  with 
the  use  of  high-speed  digital  computers  to  estimate  surface 
currents  flowing  on  the  body  by  a  finite-difference  solution 
to  a  network  of  simultaneous  equations  [Ref.  5].  Because  of 
the  limitations  of  computation  time  and  storage  capability, 
the  finite-difference  solution  is  applicable  only  when  the 
dimensions  of  the  target  do  not  exceed  a  few  wavelengths. 
For  almost  all  practical  problems,  the  use  of  approximate 
methods  is  the  only  practicable  approach. 

With  surface  current  distribution  determined,  other 
quantities  such  as  the  far  field  pattern  and  the  scattering 
cross  section  are  obtained  from  the  scattered  field  evalu- 
ated at  the  begining. 

B.   SCATTERING  BY  A  SPHERE 

A  sphere  is  the  simplest  object  susceptible  to  exact 
solution  by  the  method  of  separation  of  variables.  The 
literature  on  electromagnetic  scattering  by  spheres  is 
extensive,  and  a  great  deal  is  now  known  in  detail 
concerning   the    surface   current    distribution   and    the 
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scattered  field  in  the  neighborhood  of  the  sphere  [Ref.  6], 
The  exact  solution  of  the  scattering  of  a  plane  EM  wave  by 
an  isotropic,  homogeneous  sphere  was  obtained  by  Mie  in  1908 
and  is  usually  refered  as  the  Mie  solution. 

1.   Exact  Solution 


*.  9  . 


The  geometry  of  the  problem  in  spherical  coordinates 
is  depicted  in  the  Figure  2.2. 
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Figure  2.2    Geometry  of  Scattering  by  Sphere 

Let  us  consider  a  sphere   with  a  relative  dielectric 

constant   e    and   an   incident    wave   propagating   in   the 

r 

z-direction  and  polarized  in  x-direction. 
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In  the  spherical  coordinate  system  it  is  possible  to 
express  any  electromagnetic  field  in  terms  of  two  scalar 
functions  (Hertz  potentials)  which  are  the  radial  components 
of  the  electric  and  magnetic  Hertz  vectors  [Ref.  7], 


II  =  IIi-  r      and      II  -  n2«  r  (2.2) 

e  m 


Both  ITi  and  II 2  satisfy  a  scalar  wave  equation 


2    2 

(  V  +  k2  ) •  n  =  0       for   r  >  a  (2.3) 


Electric    and   magnetic    field    can   be    expressed    as 


E    =    VxVx(rllir)    +    i^uo  Vx  (rll2r ) 
H   =    -iueVxCrllir)    +   VxVx(rn2r) 


(2.4) 


It  can  be  shown  that  the  general  solution  of  the 
scalar  wave  equation  in  spherical  coordinates  may  be 
obtained  by  a  linear  superposition  of  the  particular  solu- 
tions ,  thus 


rll  =  r-  En   I   ITm;  (2.5) 

n=0  m=-n  n 

=   I    I      (c  ¥  (kr)  +  d  x  Ckr)}{P(m) (cosO)} 
n=o  m=-n   n  n  n  n         n 


{a  cos(m$)  +  b  sin(m$)} 
m  m 
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The  isotropic,  homogeneous  medium  will  be  considered 
to  be  a  dielectric  so  that  the  propagation  constant  k2  is 
real,  and  the  isotropic  homogeneous  sphere  is  characterized 
by  a  propagation  constant  k.  which  may  not  be  complex.  The 
ratio  of  these  quantities  defines  the  relative  refractive 
index . 


ki     miko     mi 


m  = 


C2.6) 


n\2ko     ni2 


The  incident  field,  given  in  the  equation  (2.1)  can 
be  derived  from  the  two  scalar  functions  II i  and  H2  and  are 
expressed  in  the  following  spherical  harmonic  functions: 


i           °°  Ln"1  (VI)  i 

rllx  =  £,-= ^n'x;  Y  (k2r)P1Ccos0)cos< 


(2.7) 


k*   n=1   nCn+1) 


wi     1     °°   in_1  C  2n  +  l)  i 

rllz    =  — —  I.— ^n  X)    *  (kzr)P1(cosG)sin$ 

k|   n=1   nCn+1) 


In  the  same  way,   we  write   general  expression  for  the  scat 
tered  fields  for  r>a  using  arbitrary  constants  a   and  b  : 


rn?  =— i —  ?-i C2n+1)  a  ?  (k2r )PX  (cosQ)cos$ 

k|   n=1   nCn+1) 


(2.8) 


rnf  =-- 1—  ?  -i (2n+1)  b  5  (k2r)P1(cos0)sin$ 

ki   n-X   n(n+l) 
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It  can  be  shown  that  by  applying  boundary  conditions 

we  can  get  four  linear  equations  in  the  coefficients  a  ,  b  , 

n    n 

c  ,  and  d  .    These  equations  can  be  solved  for  that  sets  of 

n  '         n  ^ 

coefficients.    Only   first  two  are   of  interest   here,   and 
these  are  given  by 


a 
n 


b   = 
n 


where 


¥  Ca)-Y  CB)  -  mf  CB)-1?  Ca) 
n      n         n      n 


5  Ca)»¥  CB)  -  m^  CB)'C  Ca) 
n      n         n      n 


m*F  Ca)  ¥  CB)  -  ¥  CB)Y  Ca) 
n      n        n     n 


mC  (a)  T  CB)  -  T  CBK  Ca) 
n      n        n     n 


C2.9) 


a    =    k2a    =    2Tra/X    =    2irmza/X 
S    =    kia    =    2iTrn.ia/Xo    =    ma 


The    scattered    fields    E      and  E       far    from    the    sphere    are    given 

0  $ 


by 


—    •    ik2r 

E0ie 

E$    = S:C0)sin$ 

k2r 


■p     .     ik2r 
E0  ie 

EQ    =      S2CQ)cos$ 

k^r 


C2.10) 


where 
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Si  CO)  =  !,  C2n+1)   {a  tt  CcosG)  +  b  t(cosQ)  } 
n=l   f     , , v    n  n  n  n 

nCn+1) 


S2(0)  =  E,  C2n+1)   {a  T  CcosQ)  +  b  v    (cosG)} 
n=l  /  , , n     n  n  n  n 


C2.ll) 


■nCn+1) 

P1 (COS0) 


d 


tt  (cosG)  =  — ,     t  (cosG)  =  — P^cosQ)     C2.12) 


n  sinQ         n         dQ  n 


2 .   Scattering  Cross  Section 

The  radar  cross  section  of  a  target  is  a  quantita- 
tive measure  of  the  ratio  of  power  density  in  the  vector 
signal  scattered  in  the  direction  of  the  receiver  to  the 
power  density  of  the  radar  wave  incident  upon  the  target 
[Ref.  8].  The  vectoral  nature  of  the  scattering  interaction 
requires  specification  of  the  transmitter  polarization  and 
receiver  polarization  with  reference  to  target  orientation 
in  three  dimensions.  Transmitter  operating  frequency  must 
also  be  specified.  The  radar  or  back  scattering  cross 
section  is  usually  defined  to  be  independent  of  range  to  the 
target,  but  it  depends  upon  target  shape  and  material,  the 
angle  at  which  the  target  is  viewed,  radar  frequency,  and 
the  polarization  of  the  transmitting  and  receiving  antennas. 
Independence  of  the  back  scattering  cross  section  upon  the 
range  to  the  target  is  valid  under  the  far-field  condition, 
i.e.,  when  the  target  is  sufficiently  far  from  the  transmit- 
ting antenna  to  justify  the  assumption  that  the  incident 
wave  is  planar  in  the  neighborhood  of  the  receiving  antenna. 
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The  conventional  definition  of   the  scattering  cross 
section  in  terms  of  the  electric  field  can  be, 


a(0,$)  =  lim   UttR2!-^.  |  (2.13) 

R-~=°        E1 


where  R  =  r  is  the  distance  from  receiving  antenna  to 
target.  Particularly  for  the  back  scattering  $  =  0  ,  0  =  II , 
and  for  the  perfectly  conducting  sphere  |  m  1  ->°0   •  Hence, 


a 
n 

= 

¥    (a) 

n 

t'  (a) 
n 

b 
n 

= 

?    (a) 

n 

C    (a) 
n 

(2.14) 


Equation  for  S2(Q)  yields 


o    (lT^    _       •    ??r    r,n-l    (n+1/2) 

S2(7r;    -    -in££-l;         . (2.15) 


5    CaH    (a) 
n  n 


because    of, 


!Cx)    =    &h   J      ,  (x)     ,  c    (x)    =    C^  H(i;(x) 

n  2  n+%  n  2  n  +  % 


TT      (-1)      =      (-l)n+1     n(!?+1)  ,  T      (-1)      =      (-l)n         n(n+1) 


n^  2  ,  unv      -,  v      ^  2 

Equation  (2.10)  now  yields, 
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ik2r 
0    "x      k2r   nSl        C  (a)c  (a) 


=  Es  =  Eo£i!H  g  (,1}n-l  Cn+1/2)   ?  (2.16) 


Finally,   back  scattering  cross  section   and  the  phase  shift 
can  be  obtained  as , 


2   *s  2         ,     E,  (-l)n'1(n+1-2) 
a   =   iim  4irR  |4_|  =  _i2L-|2zl . 1       (2.17) 

B    R-H-       Ex      k2    5  (ka)z;  (ka) 

n      n 


where  k2  =  k   ,  a  =  ka 


Z,  (-l)n"1(n+ii) 


Phase  =  -  arg{^- ; }  (2.18) 


C  CkaK  (ka) 
n      n 


For  almost  any  finite  target,  an  examination  of  the 
scattering  cross  section  as  a  function  of  frequency  leads  to 
a  separation  of  the  frequencies  into  three  regions  each  of 
which  is  characterized  by  a  particular  type  of  behavior 
[Ref.  9].  There  is,  first  of  all,  the  LF  region  where  the 
cross  section  increases  rapidly  and  smoothly  with  fourth 
power  of  frequency.  This  Reyleigh- region  behavior  is  char- 
acteristic of  any  object  whose  greatest  dimension  is  much 
smaller  than  a  wavelength.  Figure  2.3  shows  the  variation 
of  radar  cross  section  with  the  ratio  of  sphere  circumfer- 
ence to  wavelength,  2ira/X  ,  where  "a"  is  the  sphere  radius 
[Ref.  10]. 
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Increase  in  frequency,   (or  normalized  circumference 

2iTa/X,   for   sphere)   introduces   oscillations  in   the  cross 

section  curve  which   damp  rapidly  with  increasing   27ra/X  and 

2 

which  vary   about  a   mean  value  a/ira  =  1   .The   region  where 

these  oscillations  are  prominent   Cl<  2ira/X<10 )   is  called  Mie 

or   resonant   region,    as   mentioned   earlier.     Even   for 

ka  =  2Tra/X  =  43  cross  section   is  still  varying  by   ±1%  about 
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Figure  2.3    Radar  Cross  Section  of  the  Sphere 

its  mean  value. 

In  the   high-frequency  limit   (  ka>10   ),   known   as 
optical  region,   the  cross  section  settles  down  to  the  spec- 

2 

ular  value  o    -    ira   ,  the  projected  area  of  the  sphere. 


24 


III.  MEASUREMENTS  AND  RESULTS 

The  value  of  experimental  radar  cross  section  work  is 
necessary  and  quite  obvious  when  one  realizes  that  we  have 
exact  solutions  for  the  scattering  from  only  a  few  bodies, 
although  there  are,  of  course,  excellent  approximate  solu- 
tions for  several  simple  shapes.  Much  emphasis  is  placed  on 
approximate  solutions  in  studying  radar  cross  sections,  and 
experimental  results  have  been  very  important  in  helping  to 
develop  such  solutions  and  in  providing  their  validity. 

The  principal  difficulties  experienced  in  measuring 
radar  scattering  cross  sections  are  in  the  separation  of  the 
transmitted  signal  from  the  received  signal  and  the  elimina- 
tion of  unwanted  signals  that  are  reflected  from  the  back- 
ground, foreground,  and  the  target  model  support. 

Radar  cross  section  measurement  facility  at  the  Naval 
Postgraduate  School  utilizes  CW  system  that  depends  entirely 
on  cancellation  the  unwanted  signals.  The  CW  system,  prob- 
ably the  most  reliable,  is  particularly  useful  for  small 
models  that  are  measured  at  ranges  of  15  meters  or  less . 

This  chapter  deals  with  the  description  of  the  above 
system,  its  experimental  setup,  measurement  procedures,  as 
well  as  the  validation  of  the  obtained  measurement  data,  in 
great  detail. 

A.   EXPERIMENTAL  SETUP 

The  configuration  of  the  entire  system  is  shown  in 
Figure  3.1.  System  is  divided  into  three  basic  parts  as 
follows : 
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1 .   B.ange  Geometries 

The   physical  arrangement   of   the  transmitting   and 
receiving  antennas,   target  mount,   walls,   etc.   is  usually 
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Figure  3.2    Anechoic  Chamber  at  Naval  Postgraduate  School 

called  the  range  geometry.  A  variety  of  range  geometries 
have  been  developed:  these  are  characterized  by  the  way 
they  have  been  designed  to  help  eliminate  unwanted  signals 
reflected  from  the  foreground  and  background  [Ref.  11].  The 
range  geometries  commonly  used  are  free  field  ranges  with 
the  antenna  beams  horizontal,  oblique,  or  vertical;  ground 
plane  ranges,  and  image  plane  systems.  Free  space  and  image 
plane  ranges  are  commonly  used  in  both  indoor  and  outside. 

The  new  "free-field"  CW  scattering  range  at  the 
Naval  Postgraduate  School  imploys  a  metallically  shielded 
anechoic  chamber  Figure  3.2,  [Ref.  14].  It  can  also  serve 
as  the  free-field  video  impulse  scattering  range  [Ref.  13]. 
Physical  dimensions  of  the  anechoic  chamber  are  approxi- 
mately:  length  6m,  width  3m  and  height  3m.  Anechoic  chamber 
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has  the  very  important  advantage  of  being  free  from  time 
varying  interference  due  to  external  noise  such  as  atmos- 
pheric noise,  man-made  noise  (radio,  TV,  radar,  etc.),  as 
well  as  the  weather  conditions.  This  advantage  is  obtained 
at  the  expense  of  introducing  critical  problems  that  arise 
with  reflections  from  the  side  walls,  ceiling,  and  back 
wall.  These  reflections  are  reduced  by  lining  the  interior 
of  the  chamber  with  EM  waves  absorbers.  EHP-8  microwave 
absorber  made  by  Rantec  Inc. is  utilized  which  is  designed 
for  wide  range  of  frequencies  above  1.5  GHz.  The  units  are 
fabricated  of  high  quality  urethane  foam  which  is  cut  into  a 
precise  pyramidal  configuration  for  front  and  back  walls, 
and  wedge  configuration  for  lateral  walls,  floor  and 
ceiling.  Continuing  improvement  in  absorber  performance  and 
chamber  design  and  in  electronic  schemes  for  eliminating 
unwanted  reflections  have  made  it  possible  to  make  measure- 
ments on  practically  all  types  of  targets  except  those  that 
are  too  large  to  be  compatible  with  the  chamber  dimensions. 

TABLE  I 
Antenna  Specifications 

Frequency  range  4-18  GHz 

Gain  6  -  12  dB 

VSWR  (maximum)  3.2  :  1 

Isolation  <20  dB  below  5.5  GHz 

>20  db  above   5.5  GHz 

The  target  is  supported  by  a  low  cross  section  Styrofoam 
column  to  minimize  reflections  [Ref.  12].  Detailed  analysis 
of  the  support  influence  is  given  in  section  "C"  of  this 
chapter. 

The  Naval   Postgraduate  School   radar  cross   section 
range   operates   as   bistatic    system   utilizing   two   horn 
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antennas.  Both  transmitting  and  receiving  antennas  are 
mounted  on  a  removable  panel  in  the  front  wall.  Base 
aluminum  plate  that  bears  both  antennas  has  mechanism 
allowing  antenna  adjustment  in  all  three  axes.  Adjustment 
screws  are  accessible  from  outside  the  anechoic  chamber 
providing  beam  steearing  toward  target.  Antennas  are 
designed  for  both,  linear  polarizations,  horizontal  and 
vertical.  While  used  as  horizontally  polarized,  transmitter 
and  receiver,  the  other  feeds  are  terminated  with  matched 
loads  and  vice  versa  for  the  vertically  polarized  feeds. 
Basic  antenna  characteristics  are  given  in  the  Table  I. 

2 .   Instrumentation 

Equipment  used  in  the  scattering  laboratory  of  the 
Naval  Postgraduate  School,  shown  in  the  block  diagram  Figure 
3.1,  are  listed  below: 

Microcomputer  HP-85 

Synthesized  Signal  Generator  HP  8672A 

Harmonic  Frequency  Converter  HP  8411A 

Network  Analyzer  HP  8410C 

Phase-Magnitude  Display  HP  8412B 
RF  Amplifier 

DC  Power  Supply  HP  6227B 

Digital  Multimeter  HP  3466A 

Digital  Voltmeter  HP  3456A 

Digital  Voltmeter  HP  3455A 

Flexible  Disc  Drive  HP  82901M 

Plotter  HP  7225B 
HP  Interface  Bus 
Directional  Coupler 
Coaxial  Cables 

HP-85  microcomputer  provides  necessary  signal 
processing  using  developed  computer  codes,  controls  Hp  8672A 
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signal  generator,  both  HP  3456A  and  HP  3455A  digital  voltme- 
ters, as  well  as  the  HP  82901M  disc  drive  and  HP  7225B 
plotter.  Communication  between  microcomputer  and  remote 
controled  devices  takes  place  through  HP  Interface  Bus 
(HP-IB) . 

The  HP  Model  8672A  synthesized  signal  generator 
produces  RF  oscillation  in  a  frequency  range  of  2  to  18  GHz. 
The  output  is  leveled  and  calibrated  from  +3  to  -120  dBm.  AM 
and/or  FM  modulation  modes,  and  most  other  modes  can  be 
remotely  controled  using  the  HP-IB  programming  format. 

Rf  signal  from  signal  generator  enters  the  microwave 
amplifier  whose  gain  is  23  dB .  Because  of  the  low  level 
output  from  signal  generator,  amplification  is  necessary  in 
order  to  obtain  satisfactory  S/N  ratio  in  the  receiver 
output  of  the  target  return  at  higher  frequencies.  The  power 
supply  (HP  6227B)  provides  the  amplifier  with  a  DC  voltage 
and  the  digital  multimeter  (HP  3466A)  monitors  that  voltage. 

Passing  through  directional  coupler,  the  Rf  signal 
is  split  into  two.  Unattenuated  RF  portion  feeds  transmit- 
ting horn  antenna  and  the  other  portion  is  attenuated  by  43 
dB  and  feeds  the  reference  port  of  the  harmonic  frequency 
converter.  The  test  port  input  of  the  converter  is  fed  by 
the  return  signal  collected  with  the  receiving  horn  antenna. 

The  HP  8411A  harmonic  frequency  converter  and  the  HP 
8410C  network  analyzer  convert  the  RF  signal  being  measured 
to  278  kHz  signals  that  have  the  same  magnitude  and  phase 
relationships.  These  two  units  together  with  the  HP  8412B 
phase-magnitude  display,  function  as  a  phasemeter  and  a 
ratiometer  for  direct,  continuous,  simultaneous  phase  and 
magnitude  ratio  measurement  of  RF  voltages.  The  complete 
network  analyzer  measures  phase  angles  0  to  360  and  magni- 
tude ratios  in  decibels  over  a  dynamic  range  of  60  dB . 
Harmonic  frequency  converter  automatically  tracks  the 
frequency  of  the  signal  applied  to  the  reference  input. 
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The  magnitude  and  phase  DC  voltages,  taken  from  the 
phase-magnitude  display  unit,  are  measured  by  the  HP  3456A 
and  HP  3455A  digital  voltmeters  respectively.  Those  are 
versatile  digital  voltmeters  with  AC,  DC,  OHMS,  and  math 
functions,  and  offer  the  capability  of  remote  programming. 

■3 .   Targets 

Targets  used  in  this  measurement  are  several 
aluminum  and  steel  spheres  of  different  diameters  and  series 
of  right  circular  cylinders  of  various  lengths  and  diame- 
ters, made  of  thin  walled  brass  and  stainless  steel.  Spheres 
are  used  for  calibration  as  canonical  or  standard  targets  of 
known  (computed)  radar  cross  sections.  This  measurement 
reveals  correlations  of  the  theoretically  obtained  results 
using  the  Mie  series  given  by  the  equations  (2.17,  2.18). 
Cylinders,  as  primitive  models  of  the  missiles,  are  used  to 
form  necessary  data  base  for  future  work.  Adding  caps,  fins 
and  wings,  this  primitive  model  will  be  perturbed  toward  a 
real  missile  body.  Exhaustive  measurements  of  radar  cross 
section  and  phase  shift  will  be  taken  out,  and  results 
should  reveal  versatility  of  this  scheme  in  its  application 
to  target  identification  and  radar  cross  section  reduction. 

One  advantage  of  the  radar  range,  as  mentioned 
earlier,  is  the  practicability  of  testing  models  that  are 
smaller  and  cheaper  than  full-scale  targets.  For  this 
purpose  the  radar  wavelength  is  scaled  by  the  same  factor  as 
the  dimensions  of  the  model.  If  "D0"  is  any  given  dimension 
of  the  target,   and  "D  "  is   the  equivalent  dimension  of  the 

M 

model,  the  following  scaling  relation  is  employed, 


P  =  TT-=-X^  C3.1) 

M        M 
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where  X      is   the  wavelength  used  for  the   measurement  and  X . 
M  0 

is  the  wavelength  for  which  the  target  cross  section  is 
required.  At  the  same  time,  the  measured  cross  section  of 
the  model  a  should  be  multiplied  by  p2  to  give  the  target 
cross  section  a  , 


a0    =  — —   °M  -  Piou  (3.2) 


Exact  scaling  requires  the  model  conductivity  to  be  equal 
the  target  conductivity  multiplied  by  the  ratio  "p" .  This 
parameter  is  not  critical  if  the  model  conductivity  is  high 
enough . 

B.   MEASUREMENT  PROCEDURE 

Procedure  for  measuring  back  scattered  electric  fields 
from  a  target  was  suggested  in  chapter  II.  Metallic  spheres 
are  commonly  used  as  standards  since  their  cross  section  can 
be  computed,  equation  (2.17).  Prior  to  the  start  of  taking 
measurements  on  various  targets,  system  transfer  function 
has  to  be  obtained  through  a  calibration  program. 

1.   Calibration  of  the  System 

After  run  of  the  calibration  program,  microcomputer 
HP-85  remotes  and  initializes  all  programable  units,  Figure 
3.1.  Signal  generator  is  set  on  and  begins  to  generate  RF 
energy.  This  RF ,  CW  signal  undergoes  23  dB  amplification  in 
order  to  improve  S/N  ratio.  It  passes  directional  coupler 
and  is  radiated,  by  the  transmitting  antenna,  toward  the 
target  into  the  anechoic  chamber.  Signal  return  at  the  cms 
receiving  antenna,  which  is  the  vectoral  sum  of  the  target 
echo  and  the  background  radiation  is  picked  up  by  the 
receiving   antenna   which   feeds   the   "test"   port   of   the 
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harmonic  frequency  converter.  The  background  radiation 
includes  the  direct  coupling  between  the  two  antennas  and 
the  back  scattered  signal  from  the  chamber.  The  "reference" 
port  of  the  same  instrument  is  feed  with  the  RF  signal 
passed  through  another  branch  of  the  directional  coupler. 
This  signal  is  attenuated  by  43  dB  and  enters  the  frequency 
converter  at  the  about  -20dBm.  The  converter  translates  the 
RF  signals  being  measured  to  278  kHz  while  preserving  the 
same  magnitude  and  phase  relationships.  It  also  tracks  auto- 
matically any  change  in  frequency  applied  to  its  reference 
port,  and  passes  that  signal  to  the  network  analyzer.  The 
frequency  converter,  the  network  analyzer  together  with  the 
phase-magnitude  display  unit  function  as  a  phasemeter  and  a 
ratiometer  for  direct,  continuous,  simultaneous  phase  and 
magnitude  ratio  measurement  of  RF  voltages,  as  mentioned 
earlier.  Both  magnitude  and  phase  of  the  signal  return  are 
measured  by  the  HP  3456A  and  HP  3455A  digital  voltmeters 
respectively,  and  stored  in  the  HP-85  microcomputer  which 
then  computes  the  real  and  imaginary  parts  of  the  signal. 
Separation  of  the  target  echo  from  the  background  radiation 
can  be  accomplished  by  measuring  the  background  radiation 
data  as  will  be  explained  latter. 

The  procedure  repeats  "NO"  times  at  each  frequency, 
and  after  the  average  value  of  magnitude  and  phase  are 
computed,  the  signal  generator  increments  frequency  by  one 
step,  and  the  system  goes  through  new  cycle. The  block 
diagram  of  Figure  3.3  presents  various  transfer  functions  of 
the  scattering  range  [Ref.  13],  [Ref.  14].   Those  are: 

•  HT  (co  )  -  Transfer  Function  of  the  Transmitting  Circuit 

•  H  (co  )  -  Transfer  Function  of  the  Receiving  Circuit 

•  H  (co  )  -  Antenna  Coupling 

•  Hc  (co  )  -  Clutter  Signal  From  the  Anechoic  Chamber 
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•  H  (to)   -  Back  Scattered  Signal  Normalized  (  |H  (to)|2  is 
the  cross  section  and  arg[H  (w)]   is  the  phase  shift  of 

the  target) 

T 

•  H  (w)  -  Back  Scattered  Signal  of  an  Arbitrary  Target 

Besides,  X(w)  is  input  signal  provided  by  the  signal  gener- 
ator and  the  RF  amplifier.  Y(  co  )  is  reading  of  the  network 
analyzer.  The  network  analyzer  provides  the  ratio  Y(co)  'of 
the  test  channel  signal  to  the  reference  signal  which  is  a 
complex  number  at  each  frequency  step.  Note  that  every 
block,  including  X(^)  contains  noise. 
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Figure  3.3    System  Representation  with  Target  in  Place 

The  system  transfer  function  H    (cu)   and  the  target 

T  sys 

back  scattering  H  (co  )  can  be  found  by  performing  three  sepa- 
rate measurements: 

1.  without  the  target  in  place 

2.  with  the  canonical  target  in  place 

3 .  with  the  arbitrary  target  in  place 

From  the  first  measurement,  background  data  (response  from 
the  chamber  without  target  in  place)  are  obtained.  Block 
diagram  of  Figure  3.4  represents  this  measurement. 
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Figure  3.4    System  Representation  in  Background  Measurement 

Signal  return  Yi(co)  which  corresponds  to  the  scat- 
tered signal  of  the  background  can  be  expressed  by  the  equa- 
tion (3.3). 


YiCoj)    =    HCcoH    H    (u)    +    H    (u)   }H    Ceo) 

T        A         C        R 


C3.3) 


The  procedure  described  above  is  utilized  in  the 
second  measurement  as  well,  in  order  to  acquire  data  with 
the  sphere  as  canonical  target  in  place.  For  this  case 
system  response  is  labeled  as  Y2  (  bi  )  and  the  measurement  is 
represented  with  the  block  diagram  of  Figure  3.5.  Equation 
(3.4)  represents  this  configuration. 

These  two  measurements  provide  system  calibration 
because  system  transfer  function  can  be  found: 
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Y2(o))  =  H  (co)-{H  (co)  +  H^  Ceo)  +  H„Coi)}-H  (co) 

T        A         S  C        R 


C3.4) 


Subtracting  background  given  by  equation  (3.3)  from  the  data 
obtained  from  canonical  target  given  by  equation  (3.4) 
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Figure  3.5    System  Representation  with  Canonical  Target 

yields  scattered  data  from   the  sphere  (reference  standard) 
Let  it  be  Y  k  (co)  : 


Y*Cto)  =  Y2(co)  -  Y:(co)  =  H  (co)H";Cgg)H  Coj) 

T      S      R 


C3.5) 


In  equation  (3.5)   H  (co)   can  be  obtained  from  Mie  solution, 

s 

equation  (2.17).   Equation  (3.5)   yields  H    (to)   as  complex 

sy  s 
function  of  frequency  "co"  : 
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H    C„)  =  {  HT(U)  HRU)  }  "X  =  H^IV,!,)"1        (3_6) 


After  the  calibration  of  the  system  is  completed  and 

its  transfer   function  H    (to)    obtained,   we  can   find  the 

sys 

scattered  field  for  any  arbitrary  target.  To  accomplish 
this,  we  have  to  perform  a  third  measurement  with  the  target 
in  place.  System  representation  is  given  on  the  block 
diagram  of  Figure  3.3  and  its  mathematical  interpretation  is 
as  follows : 


Y3(o>)  =  HT(a>)  {HA(aO  +  HcCu)  +  H*  CuO)  HR(co)         C3.7) 


The  background  data   Y:(co)   are  subtracted    from  the  target 
data  Y  ( co)  which  we  call  Y5(co), 


Y5Cco)    =    Y3Cco)    -   YiCw)    =    H    CaOHgCaOH    (03)  (3.8) 


T 

Target  transfer  function   H  (co)is  obtained  by  simply  pluging 
equation  (3.6)  into  equation  (3.8).  It  yields  : 


H^(cq)  =  Y5(oj)  H    (w)  C3.9) 

0  sys 


Radar  or  back  scattering  cross  section  can  be  computed  as 
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2  2 

a   =  |HT(co)  |  =  |Y5(co)  H    CaO  |  (3.10) 

B      s  sys     ' 

C 
where  |H  (  co  )  |  2  =  a     is  theoretically  computed   back  scat- 

S  TH 

tering  cross  section   of  the  sphere  used   as  reference  stan- 
dard. The  phase  shift  of  the  target  is  : 


Phase  Shift  =  arg(Y5  CoO)  +  arg{H    CuO)  (3.11) 

°  sys 


2.   Measurement  on  Spheres  and  Cylinders 

Measurement  procedures  for  both  calibration  and  a 
target  data  acquisition  are  similar  as  was  described  in 
previous  section  and  consists  of  three  steps.  After  calibra- 
tion data  from  several  spheres  are  obtained  and  compared 
with  the  computed  results  of  the  respective  target.  For  that 
purpose,  theoretical  results  were  obtained  using  Mie  solu- 
tion, given  by  the  equation  (2.17),  and  represented  in 
graphical  form  ,  Figures  3.7  to  3.13.  Computer  plots  of 
target  cross  section  and  phase  were  taken  for  three  spheres 
of  different  diameters,  in  the  frequency  range  of  2  to  18 
GHz,  incrementing  frequency  in  steps  of  .02  GHz.  For  the 
first  part  of  the  experimental  work,  the  targets  used  were 
spheres  and  in  the  second  part  were  cylinders .  All  targets 
are  labeled  with  numbers  1  through  8  for  convenience. 
Besides,  brief  description  of  each  target  is  given  in  the 
header  of  the  corresponding  graph.  It  contains  basic  phys- 
ical dimensions  of  the  target,  and  corresponding  measurement 
specifications  such  as  the  polarization  of  the  antennas  and 
the  sphere  used  for  calibration.  A  brief  review  of  the 
targets  characteristics  is  given  in  Table  II.   All  graphs  of 
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theoretical  computations  and  experimental   results  are  given 
at  the  end  of  this  chapter. 


TARGET 1 
TARGET2 
TARGET3 
TARGET4 
TARGET5 
TARGET 6 
TARGET 7 
TARGET8 


TABLE  II 
Target  Characteristics 


Sphere 

Sphere 

Sphere 

Cylinder 

Cylinder 

Cylinder 

Cylinder 

Cylinder 


made  of  alum. , 
made  of  steel , 
made  of  steel , 
made  of  brass , 
made  of  steel, 
made  of  brass , 
made  of  brass , 
made  of  brass , 


D  =  .875" 
D  =  1.125' 
D  =  3.187' 


D 

= 

.75" 

L 

= 

3.875 

D 

= 

.75" 

,    L 

= 

3.875 

D 

= 

.5" 

,    L 

= 

3". 

D 

= 

.375" 

,    L 

= 

1.875 

D 

= 

.75" 

■    L 

= 

6". 

Measurements  on  spheres  and  cylinders  were  done  for 
a  frequency  range  of  10  to  15  GHz  in  steps  of  0.1  GHz. 
Graphs  of  target  cross  section  and  phase  versus  frequency 
are  represented  by  dots  overlayed  with  corresponding  theo- 
retical curves. 

C.   VALIDATION  OF  MEASUREMENT  DATA 

Recognition  of  sources  of  errors  in  measurement  data  is 
usually  difficult.  In  a  few  special  cases,  it  is  possible  to 
recognize  the  presence  of  error  and  determine  its  source  by 
observing  the  deviation  of  a  target  cross  section  vs. 
frequency  plot  from  anticipated  behavior.  The  most  common 
difficulties  are  system  noise  and  target-range  problems.  The 
latter  include  1).  the  requirement  that  the  field  incident 
on  the  target  (model)  be  a  reasonable  approximation  to  a 
uniform  plane  wave  over  the  volume   of  space  occupied  by  the 


39 


scattering  object,  2).  the  interference  of  extraneous 
reflections  from  objects  other  than  target,  and  3).  method 
of  supporting  the  target  [Ref.  12],  [Ref.  15].  Each  of 
these  problems  as  well  as  their  effect  to  the  measurement 
results  will  be  discussed  separately. 

1 .  System  Noise 

The  main  sources  of  the  system  noise  come  from  the 
RF  amplifier,  the  network  analyzer  and  from  thermal  noise  in 
the  Rf  signal  path  and  the  electronic  circuitry.  By  using 
the  anechoic  chamber  which  is  metallically  shielded,  almost 
all  external  noises  were  eliminated.  The  noise  from  the  Rf 
amplifier  is  present  in  Y(co)  because  the  test  channel  signal 
is  a  combination  of  signals  output  from  the  amplifier  at 
different  instances  while  the  reference  channel  signal  left 
the  amplifier  at  a  particular  instance.  Thus  their  ratio 
will  contain  time  varying  noises. The  network  analyzer  is  a 
noise  source  because  for  every  set  of  target  data  obtained, 
the  generator  steps  through  the  frequency  several  times. At 
one  frequency,  the  frequency  converter  may  not  pick  the  same 
harmonics  each  time  the  generator  steps  on  this 
frequency . This  causes  variations  in  the  network  analyzer 
responses.  The  other  noises  are  thermal  noises.  All  these 
components  of  the  system  noise  sources  produce  random  noises 
and  can  be  significantly  reduced  by  averaging.  Averaging 
process  was  done  twice,  first  during  the  acquisition  of  the 
magnitude  and  phase  data  at  a  single  frequency.  The  entire 
measurement  procedure  is  then  repeated  three  or  four  times 
and  averaged. 

2 .  Antenna- Target  Distance 

In  general,  the  fields  scattered  by  an  object  depend 
on  the  variation,  in  phase  and  magnitude,  of  the  incident 
field   over  the   surface  of   the  object   as  well   as  on   the 
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object's  material,  size,  and  shape  [Ref.  16].  The  scat- 
tering cross  section  of  a  target  is  defined  for  a  plane 
incident  wave  and  for  the  far  field  scattered  by  the  target. 
Because  of  the  confines  of  an  indoor  range,  non-plane  wave 
illumination  of  the  target  is  a  significant  source  of  meas- 
urement error.  The  far  from  an  antenna  is  usually  spherical 
along  its  boresight  axis  and  can  be  only  an  approximation  to 
a  plane  wave.  Such  a  spherical  wavef ront ,  when  referenced  to 
a  plane  perpendicular  to  the  antenna  boresight  and  tangent 
to  the  wavefront  at  the  nominal  target  position,  will 
produce,  to  the  lowest  order  of  approximation,  a  phase  devi- 
ating quadratically  from  that  of  an  ideal  plane  wave  along 
directions  perpendicular  to  the  central  ray . The  plane  wave 
approximation  to  the  far  field  is  satisfactorily  estab- 
lished, as  the  rule  of  thumb,  for  antenna- target  distance 
equal  to  or  greater  than  2D2 /A,  where  D  is  the  largest 
dimension  of  the  target  (model).  This  condition  as  usually 
imposed  relates  to  the  uniformity  of  the  phase  front  over 
the  lateral  region  occupied  by  the  target.  To  preserve  both 
magnitude  and  phase  uniformity  within  the  test  region,  the 
magnitude  gradient  along  the  antenna  boresight  should  be 
reduced,  as  much  as  possible.  Due  to  reciprocity,  the  same 
rule  of  thumb  of  2D2/A  as  the  range  for  far  field  applies  to 
the  scattered  field  from  the  target.  To  avoid  these  diffi- 
culties, a  more  stringent  far-field  condition  should  be 
applied.  In  other  words,  target  should  be  moved  farther  from 
the  transmitting  antenna. 

Initially  measurements  were  done  with  antenna  to 
target  range  of  1.5m,  that  caused  disagreement  with  theory 
for  the  spheres  6"  in  diameter  and  larger.  Distance  antenna- 
to-target  was  changed  to  3m,  but  the  influence  of  the  back- 
ground, that  will  be  explained  latter,  as  well  as  S/N  ratio 
caused  erroneous  results.  Finally,  an  optimal  range  of  2m 
was  established  to  satisfy  most  of,  often  contradictory 
requirements . 
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3 .   Background  Radiation 

Accurate  target  cross  section  measurements  require 
the  background  signal  to  be  held  at  very  low  level. 
Background  radiation  can  be  defined  as  returns  from  objects 
in  the  vicinity  of  the  target.  Background  effects  may  be 
classified  as  direct  backscatter  and  target  scatter.  Both 
sources  of  error  are  reduced  by  arranging  radar-absorbing 
material  so  as  to  minimize  the  undesired  background  signal 
intensity.  Residual  direct  backscatter  can  be  determined, 
as  it  was  done  in  this  experiment ,  by  measuring  the  return 
in  absence  of  the  target . Target  scatter,  or  bistatic 
coupling  effects  frequently  arise  from  the  strong  target- 
scatter  lobes  in  the  forward  direction  (away  from  the 
radar).  They  characterize  many  targets  and  are  often  two  or 
three  orders  of  magnitude  greater  than  the  direct-back  scat- 
tering which  we  measure.  Tor  instance,  a  sphere  which  diam- 
eter equals  6"  has  aforward- scatter  about  25  dB  stronger 
than  the  back  scattering  cross  section  [Ref.  17].  Thus  the 
possibility  that  the  scattered  energy  from  the  target  may 
illuminate  the  surroundings  to  return  significant  energy  to 
the  radar  cannot  be  ignored.  Illumination  of  the  background 
by  target  -  scattered  energy  may  be  followed  by  reradiation 
direct  to  the  radar  or  by  scattering  back  to  the  target  and 
then  to  the  radar. 

This  target  - scatterer  effects  were  not  measured  in 
the  experiment.  Even  direct  backscatter  can  not  be  deter- 
mined quite  accurately  because  the  introduction  of  the 
target  alters  the  radar  illumination  of  the  background  as 
well,  so  the  background  cancellation  is  no  larger  complete. 
The  resulting  error  depends  upon  the  target  cross  section, 
and  the  magnitude  of  the  background  before  and  after  the 
target  introduction.  It  is  generally  found  that  the  error  is 
small  if  the  radar- to- target  distance  is  small  compared  with 
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Figure  3.6    Measurement  Error  Due  to  Background 

both  the  radar-to-background  and  target-to-background 
distances  but  that  it  can  become  large  when  the  target 
distance  is  increased.  As  mentioned  earlier,  background 
errors  were  noticed   by  taking  measurements  at   distances  of 

1.5m  and  3m. 

Figure  3 . 6  shows  the  measurement  error  that  can 
result  from  the  background  as  a  function  of  the  ratio  of 
effective  background  cross  section  to  target  cross  section 
[Ref.  18].  To  prevent  error  greater  than  ldB ,  that  corre- 
sponds to  20%  in  cross  section,  it  is  necessary  that  the 
effective  background  cross  section  be  less  than  one 
hundredth  the  cross  section  of  the  target. 
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4 .   Target  Supports 

Target  support,  as  a  structure  employed  to  control 
and  permit  measurement  of  the  target  position  and  orienta- 
tion, is  another  important  source  of  measurement  error.  The 
magnitude  of  the  error  introduced  depends  upon  the  back 
scattering  from  the  target  support  relative  to  that  of  the 
target.  The  error  is  significant  when  the  target  cross 
section  is  small.  If  interference  with  the  support  is 
suspected,  reruns  of  the  experiment  with  different  support 
configurations  frequently  reveal  changes  in  both  cross 
section  and  phase.  In  the  case  of  the  measurements 
performed  on  spheres,  when  the  exact  solutions  and  corre- 
sponding graphs  for  the  target  cross  section  and  phase  are 
known,  deviations  from  those  known  patterns  pinpoint  to  the 
errors  in  the  measurements. 

Second  effect  that  may  occur  is  modification  of  the 
incident  field  by  the  support  device.  This  effect,  since  it 
represents  coupling  between  the  target  and  its  support,  is 
impossible  to  evaluate  quantitatively. 

In  this  work,  the  effect  of  the  back  scattering  from 
the  support,  as  well  as  the  effect  of  the  field  modifica- 
tion, were  significantly  reduced  using  tapered  Styrofoam 
supports . 

Finally,  a  brief  comment  on  measurement  results  will 
be  given,  along  with  computer  plots  of  target  cross  section 
and  phase  for  all  targets,  spheres  and  cylinders. 

Measurements  on  spheres  were  done  to  obtain  system 
calibration  and  reveal  its  validity  for  three  spheres  of 
0.875",  1.125"  and  3.187".  Besides  these  measurements  have 
shown  almost  all  previously  discussed  effects  on  the  accu- 
racy of  the  measurements  itself. 

Figures  3.13,  3.14  shows  target  cross  section  and 
phase   respectively   of   the    3.187"   sphere   in   diameter. 
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Immediately  after  calibration,  that  same  sphere  was  used  as 
the  target.  Agreement  between  measurement  and  theory  is  very 
good,  remembering  that  the  effect  of  the  system  noise 
couldn't  be  eliminated  completely  as  well  as  the  target- 
scatter  background.  It  is  observed  that  the  results  change 
in  time  because  of  the  random  nature  of  the  system  noise, 
and  thence  the   background  changes.   In  the   Figures  3.15  to 

3.18  same  calibration  was  used  and  results  are  not  as  good 
as  in  previous  measurement,  but  they  still  are  within  ±8%  of 
the  theory.  The  smaller  the  target  cross  section  the  bigger 
the  deviation.  This  is  in  accord  with  the  consideration  of 
the  background,  S/N  ratio  and  support  influence  on  measure- 
ment results . 

When  the  calibration  was  performed  using  1.125" 
sphere,  and  the  data  taken  of  that  sphere  and  0.875"  sphere 
respectively,  disagreement  between  measurement  and  theory 
became  more  pronounced.  In  this  case,  smaller  cross  section 
of  the  target  decreases  S/N  ratio,  background  and  support 
effects  are  more  enhanced,  that  all  contributes  to  the 
bigger  measurement  error.   Results   are  presented  in  Figures 

3.19  to  3.22. 

Influence  of  the  target  support  reflections  and 
coupling  are  revealed  in  Figures  3.23,  3.24  in  comparison  to 
the  Figures  3.13,  3.14.  This  time  TARGET3  (3.187"  sphere) 
was  used  for  calibration  and  to  obtain  target  data  as  well, 
but  with  the  cylindrical  support  instead  of  original  one 
that  was  tapered.  Both  cross  section  and  phase  are  changed 
slightly.  If  the  target  was  sphere  of  smaller  cross  section, 
this  effect  would  have  been  more  noticable. 

Measurements  on  cylinders  were  done  using  3.187" 
sphere  for  system  calibration  (CAL3).  Same  frequency  range  , 
as  for  the  spheres,  from  10  to  15GHz  was  chosen,  with 
frequency  increments  of  0.1  GHz.  Measurement  results  are 
plotted  with  solid  lines  to   reveal  approximate  shape  of  the 
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target  cross  section  and  phase  shift.  Results  obtained  in 
several  different  measurements,  separated  in  time,  were 
reproducible  and  thence  could  be  basis  for  future  work  on 
the  perturbation  of  a  cylinder  to  a  missile  body. 

Figures  3.27  to  3.30  show  radar  cross  section  and 
phase  for  two  cylinders  of  the  same  dimensions,  but  made  of 
two  different  materials,  brass  and  stainless  steel. 
Objective  was  to  show  influence  of  conductivity  on  the 
changes  in  target  cross  section  and  phase  shift.  This  effect 
couldn't  be  observed  because  the  both  targets  are  made  of 
materials  with  comparable  ,  very  high  conductivity. 

In  the  Figures  3.31  to  3.36  targets  were  cylinders 
of  different  diameters  and  lengths.  Corresponding  graphs 
show  that  the  cross  section  patterns  depend  upon  cylinder 
diameter  only  (number  of  resonances  is  the  same),  while  the 
target  length  along  with  its  diameter  contributes  to  the 
absolute  value  of  the  cross  section.  Phase  on  the  other  hand 
depends  only  upon  the  cylinder  diameter. 

Changed  polarization  from  horizontal  to  vertical 
gave  quite  different  results  for  TARGET4  and  TARGET6  in 
comparison  to  those  of  horizontal  polarization  obtained  for 
the  same  targets,  as  was  expected.  Figures  3.36  to  3.40 
present  these  results. 
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IV.  DESCRIPTION  OF  COMPUTER  PROGRAMS 

The  computer  programs  described  in  this  chapter  are 
written  in  BASIC,  for  frequency  domain  scattering  measure- 
ments at  the  Naval  Postgraduate  School.  They  are  run  on  the 
Hewlet  Packard  -  85  microcomputer  which,  along  with  HP 
82901M  Flexible  Disc  Drive  and  HP  7225B  Plotter,  controls 
programable  instruments  in  experimental  setup  (Figure  3.1), 
performs  necessary  computations  and  store  and  plot  results 
of  the  measurements.  Four  programs  were  written  to  complete 
scattering  experiment.  First,  "SPHERE . DRIVEO"  program 
computes  theoretical  values  for  perfectly  conducting  sphere 
which  are  latter  used  for  system  calibration  by  the  program 
"CALIB. DRIVEO" .  When  the  calibration  is  completed  and 
system  transfer  function  obtained,  program  "TARGET . DRIVEO" 
serves  for  target  data  acquisition  and  computation  of  its 
back  scattering  cross  section  and  phase  shift,  and  to  plot 
those  results  immediately  after  measurement  completion. 
Measurement  results  can  be  stored  in  the  files  on  floppy 
disc,  and  if  we  want  them  to  plot  without  taking  measurement 
again,  it  is  convenient  to  use  separate  program 
"PLOT . DRIVEO" .  Programs  are  written  to  be  interactive  and 
user  friendly.  All  that  is  necessary  to  run  any  of  those 
programs  is  to  load  them  into  main  memory  and  then  issue  the 
command  "RUN". 

A.   THEORY  OF  SPHERE 

This  program,  under  the  name  "SPHERE . DRIVEO" ,  computes 
theoretical  values  of  the  back  scattering  cross  section  and 
the  phase  shift  of  a  perfectly  conducting  sphere,  using  Mie 
solution  given  by  equations  (2.17,  2.18).  Theoretical  values 
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are  stored  in  files  of  800  records,  one  for  each  frequency, 
in  the  frequency  range  from  2.02  GHz  to  18  GHz,  at  0.02  GHz 
steps.  The  main  part  of  the  program  (lines  10  through  660) 
creates  the  file  named  "H$",  that  stores  one  value  of  the 
cross  section  and  one  value  of  the  phase,  for  each  frequency 
increment  of  0.02  GHz.  It  calls  in  SUBROUTINE  I  (lines  710 
through  990),  and  this  subroutine  by  itself  calls  in 
SUBROUTINE  II  (lines  1040  through  1530),  which  together 
carry  out  the  task  of  computation  using  previously  mentioned 
equations  (2.17,  2.18).  Program  is  listed  in  following 
pages . 
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16 

2@ 
3d 
40 
50 
60 
70 
80 
90 

1  00 
1  10 

120 
130 
140 
150 
160 
170 
180 
1  90 

200 

210 
220 
230 
240 
250 
2^0 
270 
280 
290 

300 

310 
320 
330 

340 
350 
360 
370 
380 
390 

400 

410 
420 
430 
440 
450 

4  60 

470 

480 
490 

500 


i  " SPHERE . DRIVES" 


PROGRAM  TO  COMPUTE 
VALUES  OF  BACKSCATTERING 

CROSS  SECTION  AND  PHASE 
FOR  A  PERFECTLY 
CONDUCTING  SPHERE. 


THEORETICAL  VALUES  ARE 
STORED  IN  FILES  OF  808 
RECORDS,  ONE  FOR  EACH 
FREQUENCY  FROM  2.02  GHZ 
TO  18  GHZ  AT  .02  GHZ 
STEPS . 


"THEORO.DRIVEl"  FOR  THE 
.875"  DIAMETER  SPHERE 
"THE0R1  DRIVE1"  FOR  THE 
1.125"  DIAMETER  SPHERE 
-THE0R3  DRIVE1"  FOR  THE 
3.187"  DIAMETER  SPHERE 


FILE  H$  STORES  THE 
COMPUTED  RESULT 

H$="THE0R1 . DRIVE1" 
i 

A0=1  .  125*  .  0254--2  !  SPHERE 

RADIUS  IN  METERS 
i 
! 

V9=2*PI  !  PARAMETER 
i 

i 


FREQ 
GOTO 


01=2  !  S 
0.2=18  ! 
Q4=02  • 
ON  ERROR 
PURGE  H$ 
OFF  ERROR 
CREATE  H*,800 
NEW  FILE  WITH 
!  OF  16  BYTES 
!  RECORD  STRORES  ONE 
TUDE  AND  ONE  PHASE 


TARTING  FREQ  IN  GHZ 
FINAL  FREQ  IN  GHZ 

STEP  IN  GHZ 

440 


16  !  OPEN  A 
800  RECORDS 
EACH.  EVERY 
MAGNI 
DATA. 


ASSIGN*  1  TO  H* 

DIM  B8< 144 > , B9< 1 44  > , D8< 1 44 >  , 

D9U44)  !  144>L0=INT<2*K0*A0 


510 
520 


F0  =  Q1 
FOR  1=1 


TO  806 
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530 
546 
550 
560 
570 
580 
590 

600 

610 
620 
630 
640 
650 
660 
670 
680 
690 

700 

710 
720  ' 


730 
740 

750 

7*0 
770 
780 
790 

890 

810 
R20 
830 
840 
350 

360 

870 
880 
390 
900 
910 
920 
930 
940 


DISP  "FREQ  LOOP=" , I 

p  m  =  f  0  +  Q  4 

DISP  "FREQ  <GHZ;'  =  ",F0 

K1=.3-'F0  !  WAVELENGTH 

K0  =  X9.--K1     !     Wfl'x/E    NUMBER 

ensue  700 

DISP  "E  =",E8 

DISP  "P  =n,PS 

PRINT*  1,1  ;  E9.P0 

NEXT  I 

fiSSIGN#  1  TO  * 

CLEAR 

DISP  "END  OF  COMPUTATION" 

END 


SUBROUTINE  I 


L0  = 

IF 

DIS 

CUR 

Z=K 

Gns 
E8= 
E9  = 
FOR 
L  =  L 
M3= 
M9= 
A7  = 
A8= 
L>> 
A9  = 
L  >  ) 
E3  = 
E9= 
NEX 
E8= 
E9  = 
EO  = 

P0  = 

EO  = 


IN 
L0 
P 
RE 

0* 

UE 

0 
0 

N 
0- 
D8 
BS 
<L 
A7 


FOR 


T':'2*K0-tA0  +  3) 
<145  THEN  766 
"KO*A0  TOO  LARGE 
NT  ARRAY  DIM" 
AO 

1006 


=1     TO    LO 

N+l 

CL  >'-2  +  D9<L  >-^2 

<L)  -2h-B9c:l>'^2 

+  .5>.-'M8---M9 

*<B9':L)*D9<L>-B3<L>*D3<: 


A7*  < B8  < L  >  *D9  <  L ) +  B9 < L  >  tDS  < 

A8-E3 

A9-E9 

T    N 

-E8 

-E9 

E8'x2  +  E9-"2 

ATN2'CE9,E8> 

E0/KO*2*K1     !     CROSS-SECT. 


950 
960 
970 
980 
990 

100 
101 

162 
103 
104 
105 


P9=p@_vq};  i  NT  <P0-'X9  ) 
IF    PCKPI    THEN    930 
P0=P0-V9 

P0=-P0 

RETURN 


SUBROUTINE     II 
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060 

070 

090 
100 
1  10 

120 
130 

140 

150 
166 

170 
180 
190 

200 

210 
220 
230 
240 
250 
260 
270 
280 
290 
300 

310 
320 
330 
340 
350 
360 
370 
380 

390 

400 

410 
420 
430 

440 
450 

470 
480 
490 

500 

510 
520 
530 


IF  Z>L0-1  THEN  13  70 


N2=2*Z2+L8+1 

Dl=2*N2+3 

D2  =  Dl*':2*N2  +  5> 

03  =  02*  «'.2*N2  +  7> 

04=D3*<2*N2+9> 

F 1  =  1  -22-D  1  +22^-2.-  C  2*02  >  -Z2A3 

/<6*B3> 

F2=Z*<  l.-Dl-22.^D2  +  Z2^2/(2*D3 

> -22'-3- ''-b*D4>  ) 

tt=2*Z2 

S1=F1 

F1=C2*M+1 )*F1/Z-F2 

F2  =  S1 

IF  ABS-'FIXI  .  E100  THEN  1230 

F1=F1*1  E-100 

F2=F2*1 .E-100 

Sl=Sltl  E-100 

!t  =  M-l 

IF  M+1>L0  THEN  1160 

B8<L0>=F2 

B8<L9-1>=F1 

H0=L0-2 

FOR  K=l  TO  H0 

N=L0-K-1 

B8 < N >  =  < 2*N+3 > *B8 <  N+  1  >  -'Z-B8 < 

N  +  2> 

NEXT  K 

fll  =  <:SINCZ^"Z-C0SCZ>>/B8Cl) 

FOR  K=l  TO  LO 

B8»:K>=R1*B8<K> 

NEXT  K 

GOTO  1420 

B8  < 1 >  =S I N  C  Z } 'Z-COS  < Z  > 

B8<2  >  =  <3-'ZA2- 1  >*SI  N  <  Z>  -3*C0 

S  <Z ')  --'Z 

FOR  N=3  TO  L8 

B8  <  N  >  =  <  2*N- 1  >  *B8  d  N- 1  >  -'Z-B8  ( 

N-2) 

NEXT  N 

B9<l>=-SINiZ>-C0S<Z>.'Z 

B9<2>=<1-3^ZA2>*C0S<Z>-3*SI 

N<Z>/Z 

FOR  N=3  TO  L0 

B9  (  N  )  =  <  2*N-  1  )  *B9  (.  N-  1  >  -'Z-B9  < 

N-2> 

NEXT  N 

08 <  1  >  =  <  1  - 1  'Z~2 > *S I N  f. Z  :>  +COS < 

Z  >  y  z 

D9<1  }  =  <1,'Z~2-1  }*CQS<Z>+SIN<: 

Z>  'Z 

F0R^N=2  TO  L0 

D8  <  N >  =B8  <  N - 1  >  -N*B8  (  N  )  ,'Z 

09  <  N  )=B9<N-1)  -N*B9  (.  N  )  ,'Z 

NEXT  N 

RETURN 
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B.   CALIBRATION  PROGRAM 

The  calibration  program  "CALIB . DRIVEO"  computes  system 
transfer  function  using  a  perfectly  conducting  sphere 
(canonical  target),  over  the  frequency  range  of  10-15  GHz, 
at  0.1  GHz  steps.  Computation  is  based  on  theoretical  values 
using  the  program  "SPHERE . DRIVE" ,  and  the  results  are  stored 
in  the  file  "C$". 

The  main  program  (lines  10  through  1240)  starts  with 
comments  and  necessary  declarations  of  the  variables  used 
(lines  10  through  520).  Further,  it  takes  theoretical 
values  from  H$  file  and  stores  them  into  main  memory.  For 
the  documentation  of  the  results  and  user  convenience, 
HEADER  SUBROUTINE  is  called  in  (lines  2160  through  2280). 
To  compute  system  transfer  function,  background  data  as  well 
as  the  data  acquired  from  canonical  target  must  be  avail- 
able. The  main  program  ofers  two  possibilities,  whether  to 
use  the  background  data  obtained  in  the  most  recent  calibra- 
tion, or  to  take  them  immediately  after  program  is  run.  If 
user  decides  to  take  new  background  data,  program  remotes 
all  devices  establishing  their  initial  setup  (lines  530 
through  870),  and  calls  in  BACKGROUND  DATA  COLLECTION 
SUBROUTINE.  This  subroutine  serves  to  increment  frequency  of 
the  signal  generator  (device  719)  Ml=51  times  in  steps  of  F9 
=0.1  GHz,  and  in  each  loop  to  take  averaged  data  of  the 
magnitude  and  phase,  acquired  by  the  digital  voltmeters.  It 
computes  their  real  and  imaginary  parts,  to  be  stored  in  the 
file  A$  as  well  (lines  1300  through  1490).  To  do  so,  it 
calls  in  SUBROUTINE  TO  ENTER  MAGNITUDE  AND  PHASE  DATA  from 
digital  voltmeters.  HP  3456A  digital  voltmeter  (device  720) 
reads  magnitude,  and  HP  3455A  digital  voltmeter  (device  722) 
reads  phase  (Figure  3.1).  In  each  loop  both  digital  voltm- 
eters simultaneously  take  N0=10  readings  of  the  magnitude 
and  phase,  average  them  and  pass  those  values  to  the  back- 
ground data  collection  subroutine  (lines  1550  through  1830). 
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The  main  program  proceedes  to  compute  canonical  target 
data,  after  it  is  placed  into  anechoic  chamber,  by  calling 
in  TARGET  DATA  COLLECTION  SUBROUTINE  .  This  subroutine  does 
exactly  the  same  job  the  background  data  collection  subrou- 
tine does.  It  calls  in  subroutine  to  enter  magnitude  and 
phase  data  from  digital  voltmeters.  Results  are  then  aver- 
aged, converted  to  real  and  imaginary  parts  and  passed  to 
the  main  program  (lines  1890  through  2090).  Having  both 
background  and  target  data,  main  program  computes  system 
transfer  function  (lines  1040  through  1240).  That  concludes 
calibration  procedure.   The  program  is  listed  bellow. 
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10 

20 

!  "CfiLIB. DRIVES" 

j 

30 
40 

! 

1  PROGRAM  TO  COMPUTE  SYSTEM 

56 

i  TRANSFER  FUNCTION  USING  fl 

60 

!  PERFECTLY  CONDUCTING 

70 

!  SPHERE  OVER  10-15  GHZ  AT 

80 

!  .1  GHZ  STEPS.  COMPUTATION 

90 

'  BASED  ON  THE  THEORETICAL 

100 

!  VALUES  OBTAINED  USING  THE 

ne 

!  PROGRAM  "SPHERE. DRIVEO" . 

1  20 

i 

130 

!  THE  RESULTED  SYSTEM  TRANS- 

140 

!  FER  FUNCTION  IS  STORED  AS = 

150 

i 

160 

i 

170 

!  "CALIBO.DRIVE1"  <0.875"> 

1  80 

'  -CALIE1 .DRIVE  1"  U.125"> 

1  90 

!  "CALIB3.DRIVE1"  C3.1S7"> 

200 

; 

210 

! 

220 

!  C$  IS  THE  FILE  STORING 

230 

!  THE  SYS.  TRANSF.  FUNCTION 

240 

i 

250 

'  HS     IS  THE  FILE  STORING 

260 

'  THEORETICAL  DATA 

270 

OOft 

!  OF  THE  SPHERE 
i 

•L.  O  %.' 

290 

!  S$  DESCRIBES  THE  SPHERE 

300 

i 

310 

C$="CALIB3. DRIVE1" 

320 

H$="THEQR3.DRIVE1" 

33Q 

S$="3.137  INCH  SPHERE" 

340 

fl*=HBKGRND.DRIVEl" 

350 

i 

360 

■J7Q 

X9=2*PI  !  A  PARAMETER 
i 

■3  t  f 

380 

OPTION  BASE  1 

390 

N0=1O  !  DIMENSION  OF  WO 

400 

F9=.l  !  FREQ.  STEP  IN  GHZ 

410 

Ml =51  !  Ml=(U9-L9>/F9+2 

420 

L9=10  !  LOWER  FREQ. IN  GHZ 

430 

U9=15  !  UPPER  FREQ. IN  GHZ 

440 

DIM  A (51,2)  !  BACKGROUND 

DATA 

450 

DIM  B<51,2>  !  TARGET 

460 

DIM  G3<51,2>  !  THEORY 

470 

!  CREATE  C* -Ml, 16 

480 

'  CREATE  A$,M1, 16 

490 

!  STORE  CALIBRATION  DATA 

IN  A  FILE  OF  51  RECORDS 
500  !  EACH  RECORD  CONTAINS  ONE 
MAGNITUDE  AND  ONE  PHASE 
510  !  DATA  AT  A  FREQUENCY 

53Q    ASSIGN*  1  TO  H* 
540  K0=''L9-2-F9).t5O 


550  FDR  1=1  TO  Ml 

560  KR=K8+50:*F9 

579  READ*  1..K6  ;  G3C  I  ,  1  >  ,  G3<  I  ,  2  > 

536  NEXT  I 

598  ASSIGN*  1  TO  * 

686  GOSUB  21t-0 

618  DISP  "DO  YOU  WANT  TO  USE  THE 

MOST  RECENT  BACKGROUND  DATA 
?      Y-'N  " 
620  INPUT  P* 
630  IF  P£="N"  THEN  688 
640  ASSIGN*  4  TO  AS 
650  READ*  4  ;  AC, > 
660  ASSIGN*  4  TO  * 
670  GOTO  868 
680  CLEAR 
690  REMOTE  7  !  REMOTE  ALL 

DEVICES 
700  CLEAR  7  !  CLEAR  ALL  DEVICES 
710  OUTPUT  719  .:  "P",L9,  "Z1K0L3M0 

N601" 
720  CLEAR 
730  DISP  "REMOVE  TARGET  FROM  CHA 

MBER..PUSH  'CONT'  WHEN  READY" 

740  LOCAL  7 

750  BEEP  @  BEEP 

760  PAUSE 

770  REMOTE  7 

780  CLEAR 

790  DISP  "TAKING  BACKGROUND  DATA 
■I 

800  PRINT 

810  •  PRINT  "BACKGROUND  DATA" 

820  PRINT 

830  OUTPUT  719  ; "P" ,  L9, "Zl K0L3M0 

N601  " 
840  WAIT  200 
850  GOSUB  13O0 
860  CLEAR 
870  LOCAL  7 
880  DISP  "PUT  TARGET  INTO  CHAMBE 

R..  PUSH  'CONT'  WHEN  READY" 
890  DISP  "TARGET  IS  " ,  S$ 
900  BEEP  5  BEEP 
910  PAUSE 
920  REMOTE  7 
930  ASSIGN*  3  TO  A$ 
940  PRINT*  3  .;  FK,} 
950  ASSIGN*  3  TO  * 
960  CLEAR 

970  DISP  "COMPUTING  TARGET  DATA" 
980  PRINT 

990  !  PRINT  "TARGET  DATA" 
1000  PRINT 
i018  OUTPUT  719  ;  "P" ,L9,  "Z1KHL3M 

0N6O1" 
1020  WAIT  200 
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1  830 

1  040 

050 

060 

070 
080 
090 

100 

110 

120 
130 
140 
150 

160 

170 

130 
190 

200 

210 

220 
230 
240 
258 
260 
270 
230 
290 

300 

310 
320 

330 

340 

350 

370 

380 

390 
400 
410 
420 
430 
440 
450 
460 

470 

480 


GO SUB  1890 
PRINT  "  " 

PR I  NT  " TRANS .  FUNCT I  ON " , S$ 
PRINT  "  " 
RSSIGN#  2  TO  C$ 
FOR  M=l  TO  Ml 
N1=B<M..  1  >-A<M,  1> 
N2  =  EKM,2>-FKM,2> 
X6  =  G3  <  M ,  1 > • < N 1 '  2  +  N2 A2  > 
X7  =  G3  (..  M  ,  2  }  -  ATN2  <  N2 ,  N  1  > 
X7=X7-X9*INT<:X7.''X9) 
IF  X7>PI  THEN  X7=X7-X9 
PRINT*  2,M  ;  X6,X7 
!  PRINT  USING  370  ;  M,X6,X7 
IMAGE  DD. IX, "X6=" ,SD.DDDE, 1 
X, "X7=" , SO .DDDE 
NEXT  M 

ASSIGNtf  2  TO  * 
CLEAR 

DISP  "CALIBRATION  COMPLETED 
,DRTA  STORED  IN%C* 
BEEP  e  BEEP  e  BEEP 
LOCAL  7 
END 


BACKGROUND  DATA 
COLLECTION  SUBROUTINE 


AT 


ii  p  ii 


1490 


!  OUTPUT  18  TO  15  GHZ 

O  1  GHZ  STEPS 
J=10*<L9-F9>  !  FREQ. STARTS 
AT  18  GHZ  =186*6 . 1  GHZ 
FOR  K=l  TO  Ml  !  NUMBER  OF 
FREQUENCY  STEPS 
•J=J+10*F9 
IMAGE  1A,3Z, 14A 
OUTPUT  719  USING  1368  ; 
, J, "00Z1K0L3M0N6O1 " 
!  50  MSEC  WAIT  FOR 

FREQUENCY  TO  STABILIZE 
WAIT  58 

!  TAKE  DATA  IN  FROM  722 
GOSUB  1550 
Rl =A ltCOS <P1> 
I1=R1*SIN<P1> 
A<K, 1)=R1 
A  <  K ,  2  >  =  1 1 
!  PRINT  USING 
>  A(K,2> 
NEXT  K 
OUTPUT 
0N6O1" 
RETURN 


2848 


A  <  K ,  1  > 


19 


L9, "Z1K8L3M 
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5@@  ! 

518  ! 

520  ! 

530  ! 

550  !  SUBROUTINE  TO  ENTER 

MAGNITUDE  AND  PHASE  DATA 
FROM  DIGITAL  VOLTMETER 

570  !  PREPARE  DIG.  VOLTMETERS 

580  !  TO  SEND  MAGNITUDE 

590  !  AND  PHASE  DATA 

600  ! 

610  !  NO  READINGS  TAKEN  AND 

AVERAGED  FOR  ONE  ENTRY 

620  V1=0 

630  F1=0 

640  FOR  L=l  TO  N8 

650  OUTPUT  720  ; "POFl R 1T1Z1FLOM 

0" 

66f4    WAIT  10 

670  ENTER  720  ;  V 

680  WAIT  10 

680  OUTPUT  722  i "Fl R7T1 M3A0H1 H 

700  WAIT  10 

710  ENTER  722  ;  F 

720  V1=V1+V 

730  F1=F1+F 

740  WAIT  18 

750  NEXT  L 

760  V=V1^N0 

770  F=F1/N0 

780  A1  =  10M," 

780  P1=100*F 

800  P1=DTR<P1) 

810  !  PRINT  USING  1810  ,     K,fil,P 

1 
820  IMAGE  DD,2X, "A=" ,  MD  DDDE.2X 

, °P=",SD.DDDE 
830  RETURN 
840  ! 
850  ! 
868  ! 
878  ' 
880  ! 
890  !  TARGET  DATA  COLLECTION 

SUBROUTINE 
900  ! 
910  !  OUTPUT  10  TO  15  GHZ 

AT  O  1  GHZ  STEPS 
920  J=10*<L9-F9>  !  FREQ . STARTS 

AT  10  GHZ  =100*0.1  GHZ 
930  FOP  K=l  TO  Ml  !  NUMBER  OF 
FREQUENCY  STEPS 
940  J  =  J+10:*:F9 

950  IMAGE  1A,3Z, 14A 
960  OUTPUT  719  USING  1950  ,  "P" 
■J, "00Z1K0L3M0N6O1" 
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1976 

1938 
199S 

2009 
2010 

2020 
2030 
2040 
2050 

2060 

2070 
20S0 

2090 

2100 

2110 
2120 
2130 
2140 
2150 
2160 
2170 
21  30 
2  1  90 

2200 
2210 
2220 

2230 
2240 
2250 
2266 

2270 
2230 


i  50  MSEC  WRIT  FOR 

!  FREQUENCY  TO  STABILIZE 

i  TAKE  DfiTfi  IN  FROM  722 

GOSUB  1550 

R1=A1.¥C0S<P1> 

I1=A1*SIN'CP1) 

EKK,  1>=R1 

B<K,2>=I1 

!  PRINT  USING  2840  ;  B<K,1> 

,BCK,2> 

I MflGE  4X , " R= " , SD . DDDE , 2X , " I 

=  "  ,SD.DDDE 

NEXT  K 

OUTPUT  719  ;  "P" ,L9, "21K0L3M 

0N6O1 " 

RETURN 


HERDER  SUBROUTINE 


PRINT  "  " 

PRINT  "  " 

CLERR 

DISP  "CALIBRATION  STANDARD" 

,  St 

DISP  "ENTER  TGT  RNG  IN  M" 

INPUT  Rt 

PRINT  "CRLIBRRTION  STRNDRRD 

"  ,  S* 

PRINT  "TARGET  RANGE", R* 

PRINT  "***.*#*******#*****" 

PRINT  "*:*:****************:» 

PRINT  "  " 

CLEAR 

RETURN 
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C.   TARGET  PROGRAM 

This  program  computes  back  scattering  cross  section  and 
phase  of  the  target  signal  return,  using  calibration  data 
obtained  in  the  program  "CALIB . DRIVEO"  and  stored  in  the 
file  C$ .  Results  are  stored  in  the  file  G$ .  Comments  and 
necessary  program  variable  declarations  are  given  at  the 
begining  of  the  main  program  (lines  10  through  400). 
Calibration  results  are  then  read  into  computer  memory,  as 
well  as  the  background  data  (lines  410  through  550).  If  the 
user  wants  to  take  new  background  data,  program  calls  in 
background  data  collection  subroutine  and  the  subroutine  to 
enter  magnitude  and  phase  data  from  digital  voltmeters. 
Target  is  then  placed  into  the  anechoic  chamber  and  the 
target  data  collection  subroutine  is  called  along  with 
subroutine  to  enter  magnitude  and  phase  data  from  digital 
voltmeters.  Procedures  for  both  background  and  target  data 
collection  are  described  in  the  calibration  program  section 
in  detail.  Having  calibration,  background  and  target  data 
collected,  the  main  program  computes  target  cross  section 
and  phase  shift  (lines  810  through  930).  Cross  section  and 
phase  shift  data  could  be  printed,  plotted  or  both,  stored 
in  a  newly  created  or  an  existing  file.  New  results  could  be 
averaged  with  the  measurement  results  obtained  and  stored 
earlier.  This  procedure  allows  us  to  reduce  system  noise, 
and  can  be  repeated  as  many  times  as  we  want  to  (lines  940 
through  2100).  Header  subroutine  (lines  2990  through  3190) 
is  written  to  improve  data  readability  and  to  provide  their 
documentation. 

MAGNITUDE  PLOTTING  SUBROUTINE  and  PHASE  PLOTTING 
SUBROUTINE  serve  to  plot  target  cross  section  and  phase 
shift  respectively,  from  the  results  obtained  in  this 
program,  using  HP  7225B  plotter  (lines  3230  through  4960). 
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On  the  other  hand,  if  we  want  to  plot  earlier  obtained 
and  stored  results,  "PLOT . DRIVEO"  program  has  to  be  used. 
This  program  is  developed  using  above  mentioned  plotting 
subroutines.   Target  and  plotting  program  are  listed  bellow. 
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10 

29 

3@ 
40 
59 
60 
70 
89 
90 

1  @0 
1  10 

120 
130 
140 

150 
160 

170 
186 

190 

200 

210 
220 
230 
240 
250 
260 
270 
280 
290 

300 

310 
320 
330 
340 
350 

3£fl 
370 
330 
390 

400 

410 
420 
430 
440 
450 
460 

470 
480 

490 


"TARGET .DRIVES 


PROGRAM  TO  COMPUTE  TARGET 
BACKSCATTERING  CROSS 
SECTION  AND  PHASE  USING 
C$  DATA  AND  TO  STORE 
RESULTS  IN  THE  FILE  G$ . 


FREQUENCIES: 10-15  GHZ  AT 
0  .  1  GHZ  STEPS 


FILE  C*  STORES  THE  SYSTEM 
TRANSFER  FUNCTION 
t 

!  FILE  G*  STORES  TARGET  DATA 

OBTAINED  FROM  THIS  PROGRAM 
! 

!  FILE  H*  STORES  THEORETICAL 
VALUES  FOR  PLOTTING  OVERLAY 
i 

C$="CALIB3.DRIVE1" 
G$=nSCR.DRIVEl" 
H*= "THE0R3. DRIVE  1 " 
A*=HBKGRND.DRIVE1" 

CREATE  G$,52,24 

CREATE  At ,51, 16 


OPTION 

N0=10 

DIM  A< 

DIM 
DIM 
DIM 
DIM 


BASE  1 
!  DIMENSION 
5 1 ,  2 '? 


OF  WO 


G4 
NC 
M9 


HI  =51 

F9=.  1 
U9=15 
L9=10 
DIM  T< 
TH 
X9=2*P 
ASSIGN 
READ# 
ASSIGN 
!  MAT 
IMAGE 
REMOTE 

CLEAR 
OUTPUT 
!  INI 
CLEAR 


51, 
■;"5l 
51, 
<51 
!  N 

!  F 
!  U 
!  L 

800 

EOR 
I 

#  1 
1  , 

#  1 
PRI 


2>  !  TARGET 

,2)  !  CALIBRATION 

3)  !  RESULT 

,  3) 

UMBER  OF  FREQ. 

CHECKED 
REQ . STEPS 
PPER  FREQ.  IN  GHZ 
OWER  FREQ.  IN  GHZ 
,2>  !  STORES 
ETICAL  DATA 


7  ! 

71 

TIA 


TO  C* 

G4<  ,  > 

TO  * 
NT  USING  330  .;  G4 
3D  .  4D 

!  REMOTE  ALL 
DEVICES 

CLEAR  ALL  DEVICES 
9  ;  "P1Z1K0L3M0N6O1 
L  SETUP  OF  719 
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5fifi  DISP  "DO  YOU  WANT  TO  USE  THE 
MOST  RECENT  BACKGROUND 
DflTR?     WN" 
510  INPUT  P* 
520  IF  P*="M"  THEN  576 

539  ASSIGN*  4  TO  fl* 

540  READ*  4  ;  A<, ) 
550  ASSIGN*  4  TO  * 
Fi£0  GOTO  70O 

570  DISP  "REMOVE  TGT  FROM  CHAMBE 

R  PUSH  'CONT'  WHEN  READY" 
530  LOCAL  7 
590  BEEP  5  BEEP 
680  PAUSE 

610  DISP  "TAKING  BKGRND  DATA" 
620  REMOTE  7 
630  OUTPUT  719  ; "P" > L9, "Z1K0L3M0 

N601"  !  INITIAL  SETUP  OF  719 
640  WRIT  100 
656  GOSUB  2166 
666    ASSIGN*  5  TO  fi$ 
670  PRINT*  5  ■  R<,> 
680  ASSIGN*  5  TO  * 
690  CLEAR 
700  DISP  "PUT  TARGET  INTO  CHAMBE 

R..PUSH  'CONT'  WHEN  READY" 
710  LOCAL  7 
720  BEEP  e  BEEP 
730  PAUSE 
740  REMOTE  7 
750  OUTPUT  719  .;  "P"  ,  L9,  "  71  K0L3M0 

N601"  !  INITIAL  SETUP  OF  719 
760  WAIT  500 
770  GOSUB  3030 
780  CLEAR 

790  DISP  "COMPUTING  TARGET  DATA" 
800  GOSUB  2730 
S10  F0=L9-F9 
S20  FOR  M=l  TO  Nl 
830  F0=FO+F9 
840  N<M, i>=F0 
850  X7=B<M, i>-fi<M, 1> 
860  X8=B<M,2>-A(M,2> 
870  X6=<X7-2+X8-2>*G4<M, 1> 
380  N<M,2>=X6 

890  XS  =  ATN2-rX8,X7)+G4<M,2> 
900  X3=X3-X9* I  NT  < X8'X9  > 
910  IF  X8>PI  THEN  X8=X8-X9 
920  NCM,3>=X3 
930  NEXT  M 

940  DISP  "PRINT  DATA"?  Y/H" 
950  BEEP  »?  BEEP 
960  INPUT  P$ 
970  IF  P$="N"  THEN  1616 
980  PRINT  "      FREQ       CRSEC 

PHASE" 
990  MAT  PRINT  USING  1OO0  ;  N 
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1  000 

IMAGE  2X..3D.4D 

101  0 

CLEAR 

1  020 

LOCAL  7 

1  030 

DISP  "PLOT  MAGNITUDE  FOR 

THIS  MEASURMENT?  Y/N" 

1  040 

INPUT  P* 

1  050 

IF  p$="H"  THEN  1090 

1060 

DISP  "SELECT  PEN.  PUSH 

'CONT'  WHEN  READY" 

1070 

PAUSE 

10S0 

GOSUB  3230 

1090 

CLEAR 

1  100 

DISP  "PLOT  PHASE  FOR  THIS 

MEASURMENT  ?  Y-'N" 

1  110 

BEEP  @  BEEP 

1  120 

INPUT  P* 

1  130 

IF  P$="N"  THEN  1170 

1140 

DISP  "SELECT  PEN.  PUSH 

'CONT'  WHEN  READY" 

1  150 

PAUSE 

1  160 

GOSUB  4340 

1170 

CLEAR 

1  180 

DISP  "DO  YOU  WANT  TO  MAKE 

AUER.  WITH  PREVIOUS  DATA?" 

1  190 

DISP  "?Y''N" 

1200 

BEEP  e  BEEP 

1210 

INPUT  Ft 

1220 

IF  p$="Yn  THEN  1420 

1230 

DISP  "DO  YOU  WANT  TO  STORE 

DATA  ?  Y.'N  " 

1240 

INPUT  P* 

1250 

IF  P*="N"  THEN  1880 

1260 

M0=1 

1270 

DISP  "DO  YOU  WANT  TO  STORE 

DATA  IN  FILE" 

1230 

DISP  G$ 

1290 

DISP  "?  Y.'N" 

1300 

INPUT  P* 

1310 

IF  p$="Y"  THEN  1380 

1320 

DISP  "ENTER  NAME  OF  THE  DAT 

A  FILE  TO  BE  USED  FOR  STORA 

GE" 

1330 

INPUT  G* 

1340 

DISP  "IS  THIS  AN  OLD  FILE 

TO  BE  UPDATED  ?  Y'N  " 

1350 

INPUT  P* 

1360 

IF  P$="Yn  THEN  1380 

1370 

CREATE  G$,52,24 

1380 

ASSIGN*  2  TO  G$ 

1390 

PRINT*  2  ;  M0,N<, ) 

1400 

ASSIGN*  2  TO  * 

1410 

GOTO  1380 

1420 

DISP  "ARE  THE  DATA  STORED 

IN  FILE" 

1430 

DISP  G* 

1440 

DISP  "?  Y/N" 

1450 

INPUT  P$ 
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460  IF  p$="Y"  THEN  1490 
470  DISP  "ENTER  NRME  OF  DATA 
FILE  TO  BE  USED  FOR 
THE  nVERRGE" 
480  INPUT  G$ 
490  RSSIGN#  6  TO  G$ 
500  READ*  6  .=  M0,M9(,) 
510  ASSIGN*  6  TO  * 
520  FOR  K=I  TO  HI 
530  f19 < K , 2 ) =F19 < K , 2 >*M6+N < K > 2 > 
540  M9 < K  ,  2 >  =  M9 CK ,  2 > • ( M0+ 1  > 
550  M9  <  K ,  3 >  =M9  <  K , 3 ) *MO+N  < K , 3 > 
560  M9  <  K  ,  3  >  =M9  (.  K ,  3  >  •-'  (  M0+ 1  > 
570  N<K, 2>=M9<K,2> 
i80  N<K,3>=M9<K,3> 
590  NEXT  K 
600  M0=M0+1 
610  ASSIGN*  7  TO  G* 
620  PRINT*  7  ;  MO,  NO ) 
€3Q    ASSIGN*  7  TO  * 
640  PRINT  "DATA  IS  AVERAGE  OF", 

MO, "MEASURMENTS" 
650  DISP  "PRINT  DfiTR?  Y/H" 
660  BEEP  i?  BEEP 
670  INPUT  Ft 
630  IF  P*="N"  THEN  1720 
690  PRINT  "      FREQ       CRSEC 

PHRSE" 
70O  MAT  PRINT  USING  1000  ,  N 
710  IMAGE  2X..3D.4D 
720  DISP  "PLOT  MAGNITUDE?   Y-'N" 
730  BEEP  @  BEEP 
740  INPUT  P* 
750  IF  P$="H"  THEN  1790 
760  DISP  "SELECT  PEN  FOR  MAGNIT 

UDE  PLOT.  PUSH  ' CONT '  WHEN 

RERDY.  " 
770  PRUSE 
780  GOSUB  3230 
790  CLEAR 

800  DISP  "PLOT  PHRSE?   Y-'N" 
810  BEEP  e  BEEP 
820  INPUT  Ft 
830  IF  P*="N"  THEN  1870 
840  DISP  "SELECT  PEN  RND  CHANGE 
PAPER  FOR  PHASE   PLOT. 
PUSH  'CONT'  WHEN  READY." 
850  PAUSE 

860  gosub  4340 

fi^Pi  FLE  AR 

880  DISP  "DO  YOU  WANT  TO" 

890  DISP  "OBTRIN  DRTfl" 

900  DISP  "FOR  fl  NEW  TARGET?" 

910  DISP  "  " 

920  DISP  "ENTER  Y-N" 

930  INPUT  P* 

940  IF  P*="N"  THEN  2070 
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1950 

1966 
197@ 

1  980 
1990 

2000 
2010 

2020 

2030 
2040 
2050 
2060 
2070 
20S0 
2090 

2100 

2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 


DISP  "DO  YOU  WANT  TO  USE 

THE  SANE  FILE" 
DISP  G* 
DISP  "TO  STORE  HEW  DRTfl? 

Y-'N" 
INPUT  P* 

IF  p$="Y"  THEN  2066 
DISP  "ENTER  NEW  FILE  NAME 

TO  STORE  TARGET  DATA" 
INPUT  G* 
DISP  "IS  THIS  AN  OLD  FILE 

TO  BE  UPDATED?  Y-'N" 
INPUT  P$ 

IF  p$="Yn  THEN  2066 
CREATE  G*,53,24 
GOTO  460 
CLEAR 

DISP  "END  OF  PROGRAM" 
BEEP  e  BEEP  6  BEEP 
END 


BACKGROUND  DATA 
COLLECTION  SUBROUTINE 

OUTPUT  10  TO  15  GHZ 
.l=10*<rL9-F9>  !  FREQ.  STARTS 
AT  10  GHZ  TO  BE  INCREM. 
AT  O.I  GHZ  STEPS 


2200  FOR  K=l  TO  Nl  !  NUMBER  OF 
FREQUENCY  STEPS 


2210 
2220 
2230 

2240 

2250 
2260 
2270 
2280 
2290 
23O0 
2310 
2320 
2330 
2340 
2350 

2360 
2370 
2380 
2390 


FROM 


J=J+16*F9 

IMAGE  1A,3Z,14A 

OUTPUT  719  USING  2220  ;  " 

, J, "00Z1K0L3M0N6O1 " 

!  50  MSEC  WAIT  FOR 

FREQUENCY  TO  STABILIZE 
WAIT  50 

!  TAKE  DATA  IN 
GOSUB  2420 
R1=A1*C0S<P1> 
I1=A1*SIN<P1> 
A<K, 1)=R1 
A  d  K ,  2  >  =  1 1 
!  PRINT  "11=" 
!  PRINT  "Rl=" 
NEXT  K 
OUTPUT  7 
0N6O1"  ! 
RETURN 
i 
! 
i 


P" 


A  <  K ,  2  ) 

fl(K,l) 


19  i "P" ,L9, "Z1K0L3M 
IHIT.  SETUP  OF  719 
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2  4  @  0 
2  4  i  0 
2426 


2430 
2440 

2450 

2460 
2470 
2430 
2490 
2506 
2510 
2520 
2530 
2540 
2550 
2560 
2570 
2580 
2590 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 

2770 

2730 
2790 
2800 
2810 

2820 

2330 
2840 
2850 
2860 
2870 
2838 


SU 
TU 


BR OUT  I 

DE  HMD 

DIGIT 


'  PREPARE 
TO 
RERDI 
ERflGED 


!  NO 

V1=0 

W1=0 

FOR 

OUTP 

WAIT 

ENTE 

WAIT 

OUTP 

WRIT 

ENTE 

y  i  =  u 

W1=W 
WRIT 
NEXT 
V0=V 

we=w 

Rl  =  l 
Pl  =  l 
!  PR 
P1=D 
PR 
RETU 


L=l  TO 
UT  720 

10 
R  720 

10 

UT  722 

10 
R  722 
1+V0 

1+WO 
10 

L 
l'NO 
1.--H0 

00*W0 

INT  "A 
TR-rpl> 
INT  "P 
RN 


NE  TO  ENTER  MAGNI 

PHASE  DRTR  FROM 
RL  VOLTMETERS 

DIGIT. VOLTMETERS 

SEND  DRTR 
NGS  TflKEH  AND 
FOR  ONE  ENTRY 


N0 

; "F1R1T1Z1FL0M0" 

}       V0 

; "F1R7T1M3A0H1H 
;  WO 


1  =  "  ,  A 1 
1=" ,P1 


TARGET  DRTR 
COLLECTION  SUBROUTINE 


OUTPUT  10  TO  15  GHZ  RT 
O . 1  GHZ  STEPS 
J=10*<L9-F9>  !  IHITIRL  FREQ 
UENCY  RT  10  GHZ=100*O.l  GHZ 
FOR  K=l  TO  Nl  !  FREQ . STEPS 
J=J+10*F9  !  .1  GHZ  INCREM. 
IMAGE  1A,3Z,14A 
OUTPUT  719  USING  2300  ;  "Pn 
, J, a00ZlK0L3M0N6Ol" 
!  50  MSEC  WRIT  FOR 

FREQUENCY  TO  STABILIZE 
WAIT  50 
!  TAKE  DATA 
GOSUB  2420 
R1=A1*C0SCP1) 
I1=A1*SIN<P1> 
B<K, 1 >=R1 


IN  FROM 


100 


2896 
2996 
2916 
2926 
2930 

2940 
2950 
2960 
2970 
2980 
2990 

3000 

3810 
3020 
7PI30 
3940 

3050 

3960 

3070 
3080 


B  <C  K  ,  2  >  =  1 1 

i  PRINT  "R1=",B<K,1> 

i  PRINT  "Il=" ,B<K,2> 

NEXT  K 

nUTPMT  719  ; "P" ,L9, "Z1K0L3M 

0N6O1"  !  IN  IT.  SETUP  OF  719 

RETURN 


HERDER  SUBROUTINE 


D$="MONTH-'DATE.''YEAR" 

PRINT  "  " 

PRINT  "  " 

CLEAR 

DISP  "ENTER  TODAY'S  DATE  - 

MONTH, DATE, YEAR" 

INPUT  D* 

DISP  "ENTER  TGT  DESCRIPTION 


3090 

3100 

3110 
3120 
3130 
3140 
7  15A 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 

•w,  yj  >_•  <0 

3340 


INPUT  T* 
DISP  "ENTER 


TGT  RANGE" 
R$ 
D* 

"TARGET  IS  ",T$ 
"TGT  RANGE", R* 


INPUT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
CLEAR 
RETURN 


MAGNITUDE  PLOTTING 
SUBROUTINE 


PLOTTER  IS  785 
LOCATE  32, 122,20,85 
FRAME 

SEARCH  FOR  MAX   &  M1N. 

S0=N< 1,2) 

S1=S6 

FOR  M=2  TO  Nl 

IF  S0>N<M,2>  THEN  SO=N<M, 


3350 

3360 
3370 

3330 
3390 


!  IF  S1<N<M,2>  THEN  S1=N':M, 

2) 

!  NEXT  M 

DISP  "ENTER  LOWER  VALUE  FOR 

MAGNITUDE  PLOTTING" 
BEEP  @  BEEP 
INPUT  SO 
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3460  DISP  "ENTER  UPPER  VALUE  FOR 

MAGNITUDE  PLOTTING" 
3410  INPUT  Si 
3420  L1  =  INT-'L9> 
3430  Ul=CEIl_aJSO 
3440  S3=LGT<S1) 
3450  S4=INT<S3)-1 
3460  S5=S3-S4 
3470  S5=INT<10'-S5>  +  1 
3430  S4=10^S4 
3490  L0=INT<S0-'S4> 
7500  IF  S5-L0<=14  THEN  3580 
3510  IF  35-L0>=50  THEN  3559 
3520  S5= . 5*S5 
3530  S4=S4*2 
3540  GOTO  3490 
3550  S5= . 2*S5 
3560  S4=5*S4 
3570  GOTO  3490 
3530  L0=S4*LO 
3590  U0=S5*S4 
3600  D0=U0-L0 
3610  SCALE  Ll,Ul,L0,Ue 
3620  FXD  0,4 
3630  LAXES  -1,S4,L1,L0 
3640  MOVE  L1,0 
3650  FOR  K=l  TO  Nl 
3660  M9=N<K, 1> 
3670  R9  =  NC"K,2> 
3630  MOVE  M9,R9 
3690  DRAW  M9,R9 
3700  !  GOSUB  4060 
3710  NEXT  K 
3720  M5=OJl+Ll>-'2 
3730  MOVE  M5,L0-  09*09 
3740  LORG  5  5  CSIZE  3,1,0 
3750  !  LABEL  "FREQUENCY  (GHz)" 
3760  MOVE  M5,L0-. 16*00 
3770  !  LABEL  "Fig.  .  .   TARGET2 

-  Cross  Section" 
37S0  MOVE  LI- .7,  .5*«.:L0  +  U0> 
3790  LDIR  PI^'2 
3800  !  LABEL  "CROSS  SECTION  <SQ. 

METERS)" 
3810  MOVE  M5,U0+. 095*00 
3820  LDIR  0 
7£70  CSIZE  3 ■ 1-0 
3840  !  LABEL' "TARGET:  SPHERE  D=l 

.125  inch." 
3850  MOVE  M5,U0+.06*D0 
3860  !  LABEL  "POLARI 2ATI0N  ■■     HORI 

ZONTAL" 

3870  MOVE  M5, U8+ . 925*00 

3880  !  LABEL  "MATERIAL =  ALUMINUM 
■i 

3890  PENUP 
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3960  DISP 


3916 
3926 
3930 
3940 

3950 
3960 
3970 
3980 

3990 

4000 
4010 

4020 
4830 
4040 
4050 
4060 
4070 
4  080 
4090 

4100 

4110 

4120 

4130 

4140 

4150 

4160 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 

4340 

43ciPi 

4360 

4370 

4380 

4390 

4460 

4410 

4420 


"OVERLAY  THEORETICAL 
CURVE?  Y.--N" 

P* 


H$ 


PEN  IF  DESIRED 
WHEN  READY . " 

H$ 


1  >  *  T  <  J ,  2  > 


INPUT 

IF  P*="N"  THEN  4190 

DISP  "IS  THE  THEORETICAL 

DATA  STORED  IN  THE  FILE" 

DISP  "?   Y/NH 

INPUT  P$ 

IF  P$="V°  THEN  4000 

DISP  "ENTER  NAME  OF  THE 

DATA  FILE  TO  BE  PLOTTED.1 

INPUT  H* 

BEEP  S  BEEP 

DISP  "CHANGE 

.  PUSH  ' CONT 

PAUSE 

ASSIGN#  3  TO 

•Jl  =  <L9-2>*50 

J2=<U9-2>*56 

FOR  J  =  J1  TO  J2 

READ#  3,  J  .=  T<J 

NEXT  J 

ASSIGN#  3  TO  * 

F0  =  L9 

R9=T<J1, 1) 

MOVE  F0..R9 

FOR  I=J1  TO  J2 

F0=F0+ . 02 

R9=T<I, 1> 

DRAW  F6,R9 

NEXT  I 

PENUP 

RETURN 


PLOT  CROSS 
MOVE  M9,R9 
CSI2E  2,  .5,0 
LABEL  "+" 

!  IMOVE   06025,. 00625 
!  I DRAW  -.0005,0 
RETURN 


PHASE  PLOTTING  SUBROUTINE 

PLOTTER  IS  765 

LOCATE  32, 122,26,85 

FRAME 

S4=  .5 

U0=  1 

L0  =  -1 

D0=U6-L8 
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4436 
4446 
4456 
4460 
4470 
4430 
4490 
4536 
4510 
4520 
4530 
4540 
4Fi5fi 
4560 
4570 
4580 

4590 

4  6  0  0 

4610 
4620 
4630 
4640 
4656 

4660 
4670 

4680 
4690 

4700 
4710 

4720 
4730 
4740 
4750 

4760 
4770 
4780 
4790 
4800 
4810 
4820 
4830 
4840 
4850 
4860 
4870 
4880 
4890 
4900 
4910 
4  920 


SCALE  Ll,Ul,L0,Uei 

FXD  0,3 

LRXES  -1 ,S4, LI ,L0 

MOVE  L1,0 

FOR  K=l  TO  HI 

M9=N<K,  1> 

R9=N<K,3>^PI 

MOVE  M9,R9 

DRAW  M9,R9 

!  GOSUB  4680 

NEXT  K 

MOVE  M5,L0-  .09*00 

LORG  5  e  CSIZE  3, 1 ,0 

!  LABEL  "FREQUENCY  (GHz)" 

MOVE  M5,L0-. 16*00 

!  LABEL  "Fi<3.  .  .   TARGET2 

-  Phase" 

MOVE  L 1  -  .  7  ,  . 5*  < LO+UO  > 

LOIR  PL '2 

LABEL  "PHASE  CPI>" 

MOVE  M5 , U0+ . 095*08 

LDIR  0 

CSIZE  3/1-8 

!  LflBEL  nTRRGET=  SPHERE  0=1 

.125  inch." 

MOVE  M5,U0+  06*D0 

!  LABEL  "POLARIZATION ^  HORI 

ZOHTAL" 

MO',.'E  M5  ,  U0+  .  025*06 

!  LABEL  "MATERIAL ■     ALUMINUM 
ii 

PENUP 

DISP  "OVERLAY  THEORETICAL 

CURVE?  Y'N" 
INPUT  P* 

IF  P*="N"  THEN  4866 
BEEP  G  BEEP 

DISP  "CHANGE  PEN  IF  DESIRED 
.  PUSH  'CONT'  WHEN  READY." 
PAUSE 
F0  =  L9 

R9=TCJ1, 2)  .'PI 
MOVE  F0.R9 
FOR  I=J1  TO  J2 
F0=F0+ . 02 
R9  =  T<I,2>-'PI 
DRAW  F0.R9 
NEXT  I 
PENUP 
RETURN 


PLOT  DOT 
MOVE  M9,R9 
CSIZE  2, .5,0 
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493G  LABEL  "*" 
4946  !  IMOVE  .  00825,.  006i 
4950  !  I DRAW  -.6805,0 
4  960  RETURN 
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10 

29 
30 
40 
50 
£fi 
70 
8  0 
90 

1  00 
1  10 
120 

136 

1  40 
150 
160 

170 

180 
1  90 

200 

210 
220 
230 
240 
250 
260 
270 
280 
290 

300 

310 
3  2  0 
330 
340 
350 
360 
370 

380 
390 

490 

410 
420 
430 
440 
450 


460 
470 
480 
490 

500 

510 


"PLOT .DRIVES" 

PROGRAM  TO  PLOT 
BfiCKSCfiTTERING  CROSS 
SECTION  fiND  PHASE  VS. 
FREQUENCY  USING 
STORED  DRTfl. 

G*="SCR4. DRIVE1" 
H*=HTHEOR0. DRIVE  1" 
OPTION  BASE  1 
DIM  N<51,3>  !  DRTfl  TO  BE 

PLOTTED 
Nl=51  !  NUMBER  OF 

FREQUENCES  CHECKED. 
F9=.l  !  FREQ. STEPS 
119=15  !  UPPER  FREQ.  IN  GHZ 
L9=18  !  LOWER  FREQ.  IN  GHZ 
DIM  T-:'800,2>  !  STORES 

THEORETICAL  DATA 
CLEAR 

DISP  "DO  YOU  WANT  TO  PLOT" 
DISP  "FROM  THE  FILE^ " 
DISP  G$ 

DISP  "ENTER  Y.-'N?" 
INPUT  PS 


IF  P*  = 


THEN  290 


NAME  OF 
PLOTTED 


G$ 
N( 
* 


Y/N 


DISP  "ENTER 
DISP  "TO  BE 
BEEP  @  BEEP 
INPUT  G* 
ASSIGN#  1  TO 
READS  1  ;  MO 
ASSIGN*  1  TO 
X9=2*PI 

DISP  "PRINT  DATA 
BEEP  e  BEEP 
INPUT  Ft 
IF  P*="N" 
PRINT  " 

PHASE" 
MAT  PRINT 
IMAGE  2X 
CLEAR 

DISP  "PLOT  MAGNITUDES- 
BEEP  e  BEEF- 
INPUT  FS 

IF  P$  =  ,1N"  THEN  436 
DISP  "SELECT  PEN  FOR 
DE  PLOT.  PUSH  'CONT' 
ADY  .  " 
PAUSE 
GOSUB  740 
CLEAR 

DISP  "PLOT  PHASE? 
BEEP  @  BEEF- 
INPUT  P$ 


FILE 


THEN  460 
FREQ 

USING  390 
3D.4D 


CRSEC 
N 


/N" 


MAGNITU 
WHEN  RE 


•  H" 
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520 
53  @ 


54@ 
550 
560 
570 
580 
590 

600 

610 
620 
630 
640 
650 
66& 
670 
680 
690 
7Ptfi 
710 
720 
730 
740 
750 
760 
770 
780 
790 

810 
820 
830 
840 

850 

360 
370 


8S 


80 
90 

00 


10 

20 
30 
40 
50 
60 
70 
80 
90 

M0 
01 

02 


IF  p$="N"  THEN  560 
DISP  "SELECT  PEN  AND  CHRNGE 
FOP  PHfiSE   PLOT .  PUSH 
WHEN  READY . " 


1800 

"DO  YOU  WANT  TO" 

"PLOT  DflTP" 

"FORM  nNOTHER  FILE? 


PRPER 

'  COHT 

PAUSE 

GO  SUB 

CLEAR 

DISP 

DISP 

DISP 

DISP  "  " 

DISP  "ENTER 

INPUT  Pf 

IF  pt-="Y"  THEN 

CLEAR 

DISP  "END  OF  PROGRAM" 

BEEP  £  BEEP  6  BEEP 

END 


Y/N" 
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MAGNITUDE  PLOTTING 
SUBROUTINE 


IS  70S 

2, 122,20,85 


MAX.  &  If  IN. 


PLOTTER 

LOCATE 

FRAME 

SEARCH  FOR 

S8=N< 1 , 2> 

Sl=SO 

FOR  M=2  TO  Nl 

IF  S0>N<M,2)  THEN 


S0=N<M 


IF  SKN(M;2)  THEN  S1=N<M,2 


NEXT  M 
ISP  "ENTER  LOWER  VALUE 

MAGNITUDE  PLOTTING" 
EEP  &  BEEP 
NPUT  SO 
ISP  "ENTER  UPPER  VALUE 

MAGNITUDE  PLOTTING" 
NPUT  SI 
1=INTCL9> 
1=CEIL<U9> 
3=LGT<S1 > 
4=INT<S3>-1 
5 = S  3  —  S  4 

5=1  NT < 10^55) +  1 
4=10^34 
0=INT<S0--S4> 
IF  S5-L0<=14  THEN  1080 
IF  S5-L0>=50  THEN  1050 

<--  c-  _  c  -ir  r^c 

CO— . J+JJ 


FOR 


FOR 
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1836  S4=S4#2 

1 040  GOTO  390 

1050  S5= . 2*S5 

1 060  S4=5tS4 

1070  GOTO  990 

1080  L0=S4*L0 

1090  U8=S5$S4 

1 100  D0=U0-LO 

11 16  SCALE  L1,U1,L8,U8 

1  120  FtfD  0,4 

1 1 30  LRXES  - 1 , S4 , L 1 , L6 

1 140  MOVE  LI, 8 

1150  FOR  K=l  TO  Nl 

1 160  M9=N<K, 1> 

1170  R9  =  N-:K,2> 

1  ISO  MOVE  M9, R9 

U90  DRAW  M9,R9 

1280  NEXT  K 

1210  M5=<Ul+Li::-.--2 

1 220  MOVE  M5 , LO- . 89#D8 

1230  LORG  5  (5  CSI2E  3,1,0 

1240  LABEL  "FREQUENCY  (GHz)H 

1250  MOVE  M5,L0- . 16*D0 

1260  LABEL  "Fig. 3. 13.   TARGETS  - 

CROSS  SECTION" 
1270  MOVE  L1-.7,  . 5*<L8+U8> 
1280  LOIR  PL -2 
1290  LABEL  "CROSS  SECTION  <SQ.ME 

TERS>" 
1 300  MOVE  M5, U0+ . 08*00 
1310  LOIR  O 
1320  CSIZE  3,1,8 
1330  LABEL  "TARGET:  SPHERE  <D=3. 

187in..'ALUMINUM>" 
1340  MOVE  M5,U0+  03^00 
1350  LABEL  "POLARIZATION:  HORIZO 

NTAL" 
1360  PENUP 
1370  DISP  "OVERLAY  THEORETICAL 

CURVE*  Y'N" 
1 380  ! 

1390  INPUT  P* 
1400  IF  P*="N"  THEN  166© 
1410  DISP  "IS  THE  THEORETICAL 

DATA  STORED  IN  THE  FILE",H* 
1420  DISP  "?   Y/'N" 
1430  INPUT  P* 
1440  IF  P*="YU  THEN  1478 
1450  DISP  "ENTER  NAME  OF  THE 

DATA  FILE  TO  BE  PLOTTED." 
1460  INPUT  Ht 
1470  BEEP  i?  BEEP 
1480  DISP  "CHANGE  PEN  IF  DESIRED 

.  PUSH  'CONT'  WHEN  READY." 
1490  PAUSE 

1500  ASSIGN*  3  TO  H$ 
1510  Jl=<L9-2>*50 
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526 
538 
548 
556 
=,£(* 
576 
580 
598 
680 
610 
620 
S3% 
640 
650 
f.f.% 
670 
680 
698 
788 
718 
728 
738 
748 
758 
768 
778 
788 
798 

see 

818 
828 
838 
848 
858 
868 
878 
888 
898 

see 

910 
920 
938 
948 
958 
968 
978 
988 
998 

2000 

2fllPi 
2828 
2838 
2848 
285R 
2868 


2878 


l2=<U9-2>*58 

Ffip  J  =  Ji  TO  J2 

READ*  3, J  i     T<J, 1>,T<J,2: 

NEXT  J 

ASSIGN*  3  TO  * 

F8  =  L9 

R9=TCJ1, 1 > 

MOVE  F8,R9 

FOR  I=J1+1  TO  J2 

F6=F8+ . 82 

R9  =  T<I  ..  I> 

DRAW  F8,R9 

NEXT  I 

PENUP 

RETURN 


PLOT  CROSS 
MOVE  119,  R9 
CSIZE  2,  5,6 
LABEL  "+" 

!  IMGVE  . 88825, . 66825 
•  I DRRW  -.8685,6 
RETURN 


PHRSE  PLOTTING  SUBROUTINE 

PLOTTER  IS  785 

LOCRTE  32, 122,28,85 

FRAME 

S4=  .5 

L1  =  INT''L9> 

U1=CEIL-:.U9) 

U6=l 

L6  =  -l 

D8=U0-L0 

M5=<Ul+Ll>/2 

SCRLE  L1,U1,L8,U6 

FXD  6,2 

LAKES  -1,S4,L1,L6 

MOVE  LI, 8 

FOR  K=l  TO  Nl 

M9  =  NO<,  1> 

R9=NcK,3>/PI 

MOVE  M9,R9 

M9,R9 

K 

M5,L6-. 89*08 

5  e  CSIZE  3,1,0 
"FREQUENCY  (GHz)" 


DRRW 
NEXT 
HOVE 

LORG 

LABEL 

MOVE  M5,L8-  16*D0 

LRBEL  "Fis. 3. 14.   TRRGET 

PHASE" 
MOVE  Ll-7, .5*<L8+U8> 


~7      _ 
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2  0  8  6 
2  @  9  0 
2106 
2110 
2120 
2130 

2140 
2150 

2160 
2170 

2180 
2190 
2200 
2210 

2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
235R 
2360 
2370 
2380 
2390 
2400 
2410 
2420 


LDIR  PI  -2 

LABEL  "PHRSE  CPI>" 

MOVE  M5 , U0+ . 08*D0 

LDIR  0 

CSIZE  3,1,0 

LRBEL  "TARGET'  SPHERE  <D=3. 

187  in.  -'ALUMINUM)" 

MOVE  N5 ■ U0+ . 0^*00 

LRBEL  "POLARIZATION;  HORIZO 

HTRL" 

PENUP 

DISP  "OVERLAY  THEORETICRL  C 

URVE?  Y-'H" 

INPUT  P$ 

IP  P$="N"  THEN  2320 

BEEP  6  BEEP 

DISP  "CHANGE  PEN  IP  DESIRED 

.  PUSH  'CUNT'  WHEN  READY." 

PAUSE 

F0  =  L9 

R9  =  T(  Jl,  2>-'PI 

MOVE  F0,R9 

FOR  I=J1+1  TO  J2 

F0=P0+ . 82 

R9  =  TCI  ,  2>-'PI 

DRAW  F8,R9 

NEXT  I 

PEHUP 

RETURN 


PLOT  DOT 
MOVE  M9,R9 
CSIZE  2,5,8 
LABEL  "*" 

!  IMOVE  .00025-  00025 
!  I DRAW  -.0805,0 
RETURN 


110 


V.  SUMMARY  AND  CONCLUSIONS 

Demands  on  target  cross  section  measurements  and  anal- 
ysis have  increased  considerably  in  recent  years.  There  are 
several  reasons  for  this: 

•  Great  attention  has  been  paid  to  cross  section  reduc- 
tion through  the  applications  of  radio  frequency 
absorbing  materials.  Returns  from  many  targets  of  prac- 
tical interest  are  suppressed. 

•  Introduction  of  nonmetallic  materials  into  target 
construction  has  further  complicated  target  scattering 
properties . 

•  Sophisticated  radar  systems,  on  the  other  hand, 
provided  additional  means  for  target  discrimination. 
For  example,  the  polarization  characteristics  of  a 
target  can  be  acquired  and  utilized. 

In  this  work  the  calibration  of  the  CW  scattering  range 
was  performed.  The  theoretical  solution  had  been  tested 
against  experimental  results  for  several  perfectly 
conducting  spheres  of  different  diameters.  Back  scattered 
fields  were  measured  for  three  spheres  and  five  finite 
length,  thin  walled,  right  circular  cylinders. 
Corresponding  back  scattering  cross  sections  and  return 
phase  shifts  were  obtained.  Differences  between  theory  and 
experiment  were  discussed  in  terms  of  noises  introduced  by 
such  sources  as  background  reflections,  system  noise, 
target-antenna  distance  , and  target  support  problems. 
Agreement  between  theory  and  measurement  data  was  good  if 
the  calibration  and  target  data  acquisition  were  performed 
in  a  short  time  interval.  Otherwise,  as  the  time  elapses, 
differences  between  theoretical  and   measurement  data  became 
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more  and   more  pronounced.     This  problem   requires  further 
study . 

Measurements  on  the  cylinders  showed  that  for  broad  side 
incidence,  the  number  and  location  of  the  resonance  frequen- 
cies are  strongly  dependent  on  the  diameter  of  the  cylinder. 
The  back  scattering  cross  section  depends  on  both  the 
cylinder  diameter  and  its  length.  The  phase  shift,  on  the 
other  hand,  depends  only  on  the  cylinder  diameter.  A  change 
in  polarization  provides  additional  cross  section  and  phase 
data.  Several  cylinders  (targets)  might  be  distinguished  by 
comparing  acquired  data  of  their  phases  and  cross  sections 
for  two  orthogonal  polarizations.  Future  works  will  inves- 
tigate applicability  and  reliability  of  this  approach  towad 
the  identification  of  targets  of  complex  shapes. 
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